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Abstract 

This treatise investigates holomorphic functions denned on the space of 
bicomplex numbers introduced by Segre. The theory of these functions is 
associated with Fueter's theory of regular, quaternionic functions. The alge- 
bras of quaternions and bicomplex numbers are developed by making use of 
so-called complex pairs. Special attention is paid to singular bicomplex num- 
bers that lack an inverse. The elementary bicomplex functions are defined 
and their properties studied. The derivative of a bicomplex function is defined 
as the limit of a fraction with nonsingular denominator. The existence of the 
derivative amounts to the validity of the complexified Cauchy-Riemann equa- 
tions, which characterize the holomorphic bicomplex functions. It is proved 
that such a function has derivatives of all orders. The bicomplex integral is 
defined as a line integral. The condition for path independence and the bi- 
complex generalizations of Cauchy's theorem and integral formula are given. 
Finally, the relationship between the bicomplex functions and different forms 
of the Laplace equation is considered. In particular, the four-dimensional 
Laplace equation is factorized using quaternionic differential operators. The 
outcome is new classes of bicomplex functions including Fueter's regular func- 
tions. It is shown that each class contains differentiable functions. 
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1 Introduction 



A common technique in mathematics is to define a new object as an aggregate consisting 
of previously defined objects, usually of the same type. A well-known example is to 
regard a space vector as consisting of three real numbers. The number concept is also 
extended with this technique. Thus, a rational number is viewed as a pair of integers and 
a complex number as a pair of real numbers. It is only natural to ask whether one can 
carry this idea a step further by defining numbers that are represented by pairs of complex 
numbers, or complex pairs as we shall call them. It turns out that both the quaternions 
and the bicomplex numbers can be conceived this way. The former were discovered by 
William Rowan Hamilton in 1843, the latter by Corrado Segre in 1892 @ |7). Both 
types of numbers can be seen as generalizations of the complex numbers and they lead 
to corresponding generalizations of the complex functions. 

The development of the theory of quaternionic functions was initiated by Rudolf 
Fueter in 1935 He proposed that in this theory the objects of interest ought to be 
the so-called regular functions, which are defined by means of a close analogue of the 
Cauchy-Riemann equations. More specifically, with q = x + yi + zj + uk a quaternion 
written in standard notation and if): q — > ip(q) a quaternionic function, ip is regular if it 
satisfies the partial differential equation 

(1) %L +ix W + W +kx W =Q 

ox ay oz au 

where x denotes quaternionic multiplication. In the complex plane the corresponding 
formula is 

2l +i .^L = 

dx dy 

which is an equivalent way of writing the Cauchy-Riemann equations and characterizes 
a complex holomorphic (analytic) function /: a — > /(a), where a = x + iy. 

In the other direction there is the theory of bicomplex functions. It has fairly recently 
been investigated by G.Baley Price Q, whose book also contains a brief history of the 
development of the main ideas behind the theory (including an account of Segre's work). 
The bicomplex functions of interest are the holomorphic ones, which are characterized 
by the fact that they are differentiable. They are almost isomorphic to the complex 
holomorphic functions, not surprisingly, because the operations of the bicomplex algebra 
are almost isomorphic to those of the complex algebra. The bicomplex algebra has an 
anomaly, however, which has a bearing on the associated function theory: there exist 
so-called singular numbers other than that lack an inverse. 

The purpose of this treatise is to give an overview of and further develop the theory 
of holomorphic bicomplex functions. We shall among others demonstrate that these 
functions can be divided into eight, slightly different classes. Some of the functions 
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satisfying ([j]) belong to them, hence Fueter's theory is related to our theory. 

Our purpose requires that we study the algebras of the bicomplex numbers and the 
quaternions first. Although both number types have several representations, we shall 
mostly represent them by complex pairs, such a pair being of the form (a, 6), where a 
and b are complex numbers. Complex pairs will in fact be used throughout the treatise, 
because in general they greatly facilitate calculations. 

In Chapter 3 we introduce the bicomplex functions that we regard as elementary: the 
exponential and logarithm function, the hyperbolic and trigonometric functions, poly- 
nomials and the quotient function. The properties of the exponential and logarithm 
function are specially important, since they throw light on the role played by the singu- 
lar numbers in bicomplex function theory. We shall show that the exponential function 
gets only nonsingular values and that, as a consequence, the logarithm function is de- 
fined only for nonsingular arguments. The periods that we obtain for the transcendental 
functions in this chapter are also of interest. 

In Chapter 4 we consider the differentiation of bicomplex functions. The deriva- 
tive is defined in the usual way as the limit of a fraction, but the limit must be taken 
for nonsingular values of the denominator only. Rendering the derivative as a complex 
pair yields two alternative formulas, the equivalence of which amounts to the validity 
of the complexified Cauchy-Riemann equations. Fulfillment of these equations is the 
main characteristic of the holomorphic bicomplex functions. We prove that such a func- 
tion possesses derivatives of all orders. In addition, we give R 4 -representations of the 
bicomplex derivative. 

In Chapter 5 we consider the integration of bicomplex functions. The basic bicomplex 
integral is defined as a line integral, which is path independent if and only if the integrand 
is the derivative of a differentiable bicomplex function. We show that a holomorphic 
bicomplex function can be expanded into a Taylor series. In order to prove the bicomplex 
versions of Cauchy's theorem and integral formula we generalize Green's theorem and 
derive the so-called twining number, the bicomplex analogue of the complex winding 
number. 

In the last chapter we clarify the relationship between bicomplex functions and the 
Laplace equation. For this purpose we first study the derivatives of complex holomorphic 
and conjugate holomorphic functions. We factorize the four-dimensional Laplace equa- 
tion using quaternionic differential operators and concepts obtained from the similar but 
simpler factorization of the two-dimensional version of the equation. This leads to the 
aforementioned classes of bicomplex functions, each class containing also differentiable 
functions. In this context we shall have the occasion to define the derivative related 
to Fueter's regular functions, a topic that has received attention in the literature. Our 
proposal is to define the derivative as the limit of a bicomplex fraction where both the 
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numerator and the denominator are regular in Fueter's sense. 

Theorems, lemmas, definitions and formulas are numbered separately within each 
chapter. The number of a theorem, lemma or definition is always preceded by the chapter 
number. Numbered formulas are referred to by the assigned number when needed in the 
same chapter, otherwise the chapter number is added. 
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2 Quaternionic and bicomplex algebra 



In this chapter we shall make extensive use of complex pairs to develop the algebras 
of quaternions and bicomplex numbers. Although both algebras generalize the complex 
algebra, we can of course not expect them to retain all the properties of this algebra: in 
the quaternionic algebra the commutativity of the multiplication operation is lost, and 
in the bicomplex algebra we must accept the existence of the singular numbers. The 
bicomplex numbers are more important for our work, but the quaternions are not unim- 
portant. They are needed in the factorization of the four-dimensional Laplace equation 
and, above all, provide us with a norm that directly can be taken over to the bicomplex 
algebra. Thus, to begin with we shall investigate the quaternions. 

2.1 Representation of quaternions 

A quaternion is in general written as a four-dimensional vector of the form 

x + yi + zj + uk 

where x,y,z,u are real numbers and i,j, k unit vectors (generalized imaginary units). 
Given two quaternions 

q = x + yi + zj + uk 
r = X + Yi + Zj + Uk 

of this kind we may add them by applying component-wise addition, formally: 

(1) q + r = {x + X) + (y + Y)i + {z + Z)j + (u + U)k 

To denote quaternionic multiplication we shall consistently use the symbol x . The well- 
known multiplication table for the unit vectors is 

(2a) ixi = -l j x j = -1 k x k = -1 

(2b) i x j = k j x k = i k x i = j 

(2c) k x j = -i j x i = -k i x k = -j 

which yields the following expression for the non-commutative product q x r: 

(3) qxr = (xX - yY - zZ - uU) + (xY + yX + zU - uZ)i + 

(xZ -yU + zX + uY)j + (xU + yZ - zY + uX)k 

We now wish to represent the quaternions q and r as pairs of complex numbers: 
q = (a, b) , a = x + yi and b = z + ui 

r = (c,d) , c = X + Yi and d = Z + Ui 
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How should the multiplication x be expressed in terms of the complex numbers a, b, c, d? 
We find 

(4) q x r = (a • c - b ■ d* , b ■ c* + a ■ d) 

where the operators • and * stand for complex multiplication and conjugation. Writing 
the right-hand side of this formula in terms of x, y, z, u and X, Y, Z, U yields 

(5) qxr = (xX - yY - zZ - uU + (xY + yX + zU - uZ)i , 

xZ - yU + zX + uY+ (xU + yZ - zY + uX)i ) 

which is another way of rendering (|^). 

It is sometimes convenient to represent the quaternion x + yi + zj + uk as a quadruple 

(x,y,z,u) 

With q = (x, y, z, u) and r — (X, Y, Z, U), the addition and multiplication rules then look 
like 

q + r = (x + X,y + Y,z + Z,u + U) 

qxr = (xX -yY - zZ -uU , xY + yX + zU - uZ , 
xZ -yU + zX + uY , xU + yZ - zY + uX) 

In the context of quaternionic algebra the four-dimensional real space R is usually called 
the quaternionic or Hamiltonian space H. It can also be regarded as the Cartesian product 
of two complex planes, or C 2 . The different views of the space R 4 enable us to classify 
the representations we have introduced for a quaternion q = (a, b) = {x + yi, z + ui): 

H-representation: q 

C 2 -representation, complex pair: (a, b) 

(6) R 4 -representation, complex pair: (x + yi, z + ui) 
R 4 -representation, vector-form: x + yi + zj + uk 
R 4 -representation, quadruple-form: (x,y,z,u) 

As a rule the C 2 -representation leads to the shortest calculations. 

2.2 Basic quaternionic algebra with complex pairs 

Let (a, b) and (c, d) be two complex pairs representing quaternions. The pairs are equal 
if and only if their components are equal, formally: 

(7) (a, b) = (c, d) = (a = c) A (6 = d) 
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Remark. The logical operators are designated by = (equivalence), =>• (implication), V 
(disjunction), A (conjunction), and -i (negation). Their priority order is from weaker to 
stronger: =, =>, V and A, V and A have the same priority. □ 

For quaternionic addition and multiplication we found in the previous section 

(8) (a, b) + (c, d) = {a + c,b + d) 

(9) (a,b)x(c,d) = (ac - bd* , be* + ad) 

With these C 2 -formulas it is not difficult to verify the following basic properties of + and 
x (using one of the R 4 -representations for the quaternions the proofs of the first three 
properties are rather laborious): 

[(a, b) x (c, d)] x (e,/) = (a, b) x [(c, d) x (e, /)] {x is associative} 

(a, b) x [(c,d) + (e, /)] = (a, b) x (c, d) + (a, b) x (e, /) {left-distribution} 
[(c,d) + (e,/)] x (a, 6) = (c, rf) x (a, 6) + (e, /) x (a, 6) {right-distribution} 

(a, 6) x (1,0) = (a, 6) {identity element of x} 

(a, 6) + (0, 0) = (a, 6) {identity element of +} 

(a, b) x (0,0) = (0,0) {zero element of x} 

The complex pairs (1,0) and (0,0) are examples of real or scalar quaternions, whose 
general form is (A, 0), where A € R. Multiplication of a quaternion by a real quaternion 
is in a sense a distributive operation, since for any (a, 6) € C 2 

(10) (A,0) x (a, 6) = (Aa, A6) 
Moreover, it is also a commutative operation or 

(A,0) x (a, 6) = (a, 6) x (A,0) 

The quaternion (A, 0) is designated by just A, because it has all the properties of the 
real number A. The x -operator can be omitted from a quaternionic product if one of the 
factors is a real quaternion. For example, the previous formula is written A(a, b) = (a, b)X. 

2.3 Quaternionic conjugation, square norm, and division 

In this section we shall introduce some concepts based on the quaternionic conjugate. 
The terminology and notation are partly taken from B , [Ol and p3| . We first recall the 
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properties of the complex conjugate. With a and b two complex numbers the conjugation 



operator * satisfies: 

(11) (a + b)* = a* +b* {* distributes over +} 

(12) (a • b)* = a* ■ b* {* distributes over •} 

(13) (a*)* = a {* is an involution} 

We shall use the symbol * for quaternionic conjugation, too. Its definition in C 2 is 

(14) (a,b)* = (a*,-b) 



The equivalent Reformulations are: 

(x + yi, z + ui)* = (x — yi,—z — ui) 
(x + yi + zj + uk)* = x — yi — zj — uk 
(x,y,z,u)* = (x,-y,-z,-u) 

Quaternionic conjugation is an involution and distributes over addition, but not over 
multiplication. Hence, of the formulas (|Tl|), (|l2|), and jl^ ) only the first and last may be 
taken over to quaternionic space: 

(15) (q + r)* = q*+r* q,reR 

(16) (q*)* = q 

The square norm N(q) of the quaternion q is defined by 

(17) N(q) = q x q* {H-representation} 

With q = (a, b) we obtain the equivalent complex pair formula, which reveals that the 
square norm is a real number: 

(18) N((a, b)) = aa* + bb* {(^-representation} 
In R 4 the same formula looks like 

(19) N((x + yi, z + ui)) = x 2 + y 2 + z 2 + u 2 {R 4 -representation} 

If we let N(a) stand for the square norm of a complex number a, too, we have N(a) = aa*, 
and hence, by ( |l8| ) 

(20) N((a, b)) = N(a) + N(b) 
For two quaternions q and r we also have 

N(q x r) = N(q)N(r) 
7 



The function N enables us to introduce further quaternionic concepts. We denote the 
inverse of a nonzero quaternion q by q f — 1 (the standard notation q~ 1 being reserved 
for another concept), and define 

(21) ? T-i = ^<f ,«*o 

With the inverse at our disposal we can introduce a division operation for two quaternions 
r and q, but due to the non-commutativity of x we must distinquish between left-division, 
(gt — 1) x r , and right-division, r x (ql — 1). 

Finally, we define the notion of the absolute value or modulus \\q\\ of q by 

(22) NI = V^(«) ,geH 

The || 1 1 -function plays an important role in the bicomplex algebra, too, as we shall soon 
see. 

2.4 Bicomplex algebra 

The space of bicomplex numbers, which we denote by B, has many properties in common 
with the Hamiltonian space. A bicomplex number q is foremost viewed as a complex pair 
(a, b), where a, b 6 C, but all the other quaternion representations (Q) apply to it, too: 

B-representation: q 

C 2 -representation, complex pair: (a, b) 

(23) R -representation, complex pair: (x + yi, z + ui) 
R 4 -representation, vector-form: x + yi + zj + uk 
R 4 -representation, quadruple-form: (x,y,z,u) 

The equality condition of two bicomplex numbers is given by (Q) and they are added 
according to (|§|). The multiplication of bicomplex numbers, however, differs from the 
quaternionic operation. If we write a complex number x + iy as a pair of real numbers 
(x,y), ordinary complex multiplication of a = (x,y) by b = (z,u) takes the form 

a ■ b = (xz — yu, yz + xu) 

Replacement of the factors a and b in this scheme by bicomplex numbers yields bicomplex 
multiplication, designated by 0. Thus, for q — (a, b) and r = (c, d) we have the definition 

(24) qQr = (a-c-b-d,b-c + a-d) 

Evaluating this formula for q — (x + iy,z + iu) and r = (X + iY, Z + ill), we obtain the 
equivalent R 4 -formula 

qQr = (xX —yY - zZ + ull + i(xY + yX - zll - uZ) , 
xZ -yU + zX ~uY + i(xU + yZ + zY + uX)) 
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The vector-form q = x + yi + zj + uk of a bicomplex number is occasionally useful. The 
unit vectors have the complex pair representations 



(25) i=(t,Q) , j = (0, 1) , fc = (0,t) 

and are multiplied according to: 

i©i=-l j'0j = -l fc0fc = l 

i0 j = fe j Q k = —i kQ i = —j 

j i = fe k Q j = —i i fc = — j 

A bicomplex scalar (A, 0), where A £ R, will be abbreviated by just A in the same way 
as a real quaternion. The identity element (1,0) and the zero element (0,0) of are 
therefore written 1 and 0. We also abbreviate the expression A (a, b) by A(a, 6) and 
note that 

(26) \{a,b) = (\a, Xb) 

The bicomplex number (c, 0), where c € C, has all the properties of the complex number 
c, yet we shall not use any abbreviation in this case. 

2.5 Further bicomplex operations 

The 0-operator is associative, distributes over + and, most significantly, is commutative. 
It therefore permits us to generalize some of the algebraic operations of the complex 
numbers. The exponentiation of a bicomplex number is performed in the customary 
fashion: 

f ?° = 1 

| q n = qQ q n ~ x if n is an integer > 1 

Bicomplex conjugation is meaningful, too. We designate it by a postfix u and give it the 
C 2 -definition: 

(28) (a,br = (a,-b) 

It is isomorphic to the definition of complex conjugation * expressed by real pairs or 
(x,y)* — (x,—y). As a result, the u -operator inherits the properties Ji~l|) , (|l2|) , and ( |l3| ) 
of*: 

(29) (q + r) u = q u + r u q,reB 

(30) ( 9 0r) u = q u Gr u 

(31) (q u r = q 

Next we define the complex square norm CN(q) of the bicomplex number q: 

(32) CN(q) =q(Dq u {B-representation} 



(27) 



9 



a formula of the same type as (|17j). With q = (a, b) we get the C 2 -representation: 

(33) CN[{a,b)) = (a 2 + b 2 ,0) 

It shows that CN(q) is not a real number as N(q) but a complex number (which we 
nonetheless prefer to write in bicomplex notation). C/V((a,6)) is zero if and only if 
a 2 + b 2 — 0, which is the case for b = ia or b — —ia. Viewed as vectors a and b are then 
of equal length and perpendicular to each other. We thus have the relation: 

(34) CN((a, b)) = = (b = ia) V (b = -ia) 

The condition CN(q) = is crucially important in bicomplex algebra and has no essential 
analogue in complex algebra. We shall use the following terminology : 

Definition 2.1 A bicomplex number q is called singular if CN(q) = 0, otherwise it is 
called nonsingular. □ 

From the point of view of developing a bicomplex function theory it is comforting to 
observe that the bicomplex space contains simply connected domains that are not single- 
ton sets and contain merely nonsingular points. Such domains are consequently termed 
nonsingular domains. 

Example. The domain that consists of all bicomplex numbers (x + iy, z + iu), such that 
ie>0A?/>0Az>0Au>0 holds, is nonsingular. □ 

For two bicomplex numbers q and r the complex square norm satisfies 

(35) CN(q r) = CN(q) CN(r) 

from which it follows that the bicomplex product preserves nonsingularity if and only if 
both factors are nonsingular. 

The singular bicomplex numbers (a, ia) and (a, — ia) are so-called zero divisors, be- 
cause (a, ia) (a, — ia) = p5| . Hence, although g = V r = => qQr = holds for 
all q and r, the implication in the other direction is not generally valid. 

As the inverse of CN((a, b)) we can take 

( 36 ) 7^777 ,» = 1 9 1 ,9 > ) ,a 2 + b 2 ^Q 
y ' CN((a,b)) \a 2 + b 2 J r 

The inverse q^ 1 of the bicomplex number q is now defined by 

(37) ^ = cW) GqU 'M®* 
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The definition is of the same form as (|2l]). Note that only a nonsingular bicomplex 
number can possess an inverse. Inserting q = (a, b) in p^ ) yields upon application of 
(EJ), (Eq) and (pq) the equivalent C 2 -representation 



At this point the definition of the division of two bicomplex numbers r and q is straight- 
forward: 

(39) -=rQq- 1 ,CN(q)^0 

q 

The notation ^ is meaningful due to the commutativity of 0. In C 2 the same formula 
looks like 

(c, d) I c ■ a + d ■ b d ■ a — c • b 



(40) i^n^T^^T^j ' + ^° 

Using the inverse (^) we are also able to extend the definition of exponentiation (|27j) to 
negative integer exponents by 

(41) q- n = Or 1 )" , n > and CN(q) ^ 

The definitions of the inverse and of the division operation are thus based on the complex 
square norm CN . The definition of the absolute value (modulus) of a bicomplex number, 
however, cannot be based on CN; for this purpose we have to borrow formula ( ^2|) from 
quaternionic algebra: 

(42) h\\ = y/Nti) ,g£B 



Due to ( pL9| ) the |j ||-function may be identified with the Euclidean norm, which is known 
to have the following properties: 

\\q\\ > for all q £ B 

||?||=0 = ? = 

||Ag|| =| A Hkll for all A G R and q eB 

h + r\\ < \\q\\ + \\r\\ for all q,r S B 
Moreover, formula (|2^) gives us for all bicomplex (a, 6) 

II («, 6)11 < 1*1 + 1*1 

The notation | c | stands for the absolute value \J N(c) of the complex number c. 

If we visualize the bicomplex number q as a vector in R 4 , the absolute value ||g|| is 
the length of the vector. 
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2.6 Octonions and tricomplex numbers 

It remains to say something about the generalization of the quaternions and the bicom- 
plex numbers in the form of octonions and tricomplex numbers, respectively. 

An octonion can be expressed as a pair of quaternions {q, r), thus making the octo- 
nionic space O 8-dimensional. Octonionic addition is performed in the normal way by 
component-wise addition, i.e. (q, r) + (s, t) = (q + s,r + 1). Octonionic multiplication 
is defined in terms of quaternionic multiplication by 

(43a) (q, r) (s, t) = (q x s - t* x r , r x s* + t x q) q,r,s,t 6 H 

This formula should be compared with the formulas we have employed for complex and 
quaternionic multiplication: 

(43b) (x, y) ■ (z, u) — (xz — uy, yz + ux) x,y,z,u € R 

(43c) (a, b) X (c, d) = (a ■ c - d* ■ b , b ■ c* + d ■ a) a, b, c, d e C 

If the conjugate of a real number is taken to be the number itself, (43b) could be written 
(x, y) ■ (z, u) = (xz — u*y, yz* + ux). The formulas (43a), (43b) and (43c) clearly follow 
the same scheme and therefore complex numbers, quaternions and octonions have many 
properties in common. However, generalizing from quaternions to octonions causes the 
loss of the associativity of the ^-operation. Instead <E) satisfies the weaker property of 
alternation: (to (g) m) <g) n = to ® (to ® n) and to ® (n ® n) = (m ® n) (g) n for to, n G O. 

A tricomplex number, in turn, is expressed as a pair of bicomplex numbers. Addition 
of two tricomplex numbers (q, r) and (s, t) is performed by component-wise addition and 
multiplication by 

(q, r) (s, t) = (q © s - t © r , r s + 1 g) <?, r, s, < e B 

Here we have taken the liberty of overloading the symbol by applying it to the mul- 
tiplication of both bicomplex and tricomplex numbers. This operation is isomorphic to 
complex multiplication (43b), as is bicomplex multiplication (^J). 
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3 Elementary bicomplex functions 



In this chapter we study functions of the type B — > B, or bicomplex functions of a 
bicomplex variable. The benefits of using complex pairs will now become so obvious that 
the alternative to base one's reasoning entirely on one of the R 4 -representations of the 
bicomplex numbers seems exceedingly unattractive. 

We define the elementary bicomplex functions, or the exponential function, the hyper- 
bolic and trigonometric functions, polynomials, the quotient function and the logarithm 
function. These functions contain their complex counterparts as special cases, as must 
be required of them if the bicomplex function theory is to be a proper generalization of 
the corresponding theory in the complex domain H ^|, |ll|, |24| . 

3.1 The form of bicomplex functions 

Let the bicomplex variable p be given by: 

(1) P = (a,b) 

(a, b) = (x + iy, z + iu) 

A bicomplex function ip of p then has the form: 

(2) V(P) = (Ma,b),fa(a,b)) 

(<t>i(a,b),<fa(a,b)) = (^x(x,y,z,u) + iip 2 {x,y,z,u), 
ip3(x, y, z, u) + itpi(x, y, z, u)) 

The functions <f>k,k = 1,2, are of type C 2 — > C, the functions ipk,k = 1,2,3,4, of type 
R 4 -> R. 

Formulas (0) and (§) provide us with B-, C 2 - and R 4 -representations of p and ip. 
Because of our frequent reliance on the C 2 -representation it would not be inappropriate 
to look at ip as a function of two complex variables a and b rather than as a function of 
one bicomplex variable p. In any case, what is essential in the sequel is the requirement 
that the two- argument, complex functions (f>k should be holomorphic in both a and b, 
i.e. the partial derivatives and should exist. Roughly speaking, the 0fc's meet 
this condition if their values are obtained by algebraic operations and applications of 
holomorphic functions so that a and b enter in them only as complex aggregates without 
separation of their real and imaginary parts po[ . 

Example. For a = x + iy and b — z + iu the functions a 2 + b 2 and a -sin 6 are holomorphic 
in a and b, not so x 2 + y 2 + b 2 and a ■ sin(z + u). □ 
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At this point we recall some basic results from complex analysis. In general, if a 
complex function 

u(a,b) =£ 1 (x,y,z,u) + i£ 2 (x,y,z,u) 

is holomorphic in both a = x + iy and b = z + iu the partial derivative ^ is given by 
either of the formulas 

(3) ^ = 9 x £i + id x & 

oa 

(4) ^ = dyfr-idyh 

Their equivalence amounts to the truth of the Cauchy-Riemann equations 

(5) 9x6 = dy& , dx& = -d y ^ 

By making the replacements a — >■ b , x z , y — > it in (Q), (Q) and (||) we obtain the 
corresponding formulas for ^ . 

Remark. Above we abbreviated the operators J| and by d x and d y . We employ this 
notation when taking the derivative with respect to a real variable, but not when taking 
it with respect to a complex variable. 

We further left out the arguments of the functions. We allow ourselves such an 
omission whenever it is clear from the context which the arguments are. □ 

Returning to the function if) in it should be emphasized that the holomorphism of 
its component functions <f>k in a and b does not ensure that -0 itself is holomorphic in 
p. For ip to be holomorphic the ^-functions must fulfill a stronger condition, namely 
the bicomplex version of the Cauchy-Riemann equations. The issue will be discussed in 
Chapter 4. But the holomorphism of the component functions does suffice to establish 
the continuity of ip. To do so we need the limit concept lim p _» Po tp(p) — wq, which in the 
ordinary fashion is defined to be equivalent to 

(6) (Vp,e:e>0: (36 : 5 >0 : \\p-po\\ <5 => U(p) - w Q \\ < e)) 

where pgB and e,S E R. The absolute values should be computed using (2.^2|). Geo- 
metrically interpreted, the limit is taken for all p in a (deleted) neighbourhood of the 
bicomplex space or, more precisely, a four-dimensional ball of radius 6 centered at, but 
not including, the point pq. 

Falling back on (^), one can define the continuity of "0 at a point po by 

(7) lim 0(p o + Ap) = if/fa) 

Ap^O 
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The holomorphism of the functions </>&, k = 1,2, at the point (ao,&o) means that they 

d<Pk artr \ d±k 

-g^- ana 



possess partial derivatives and ^fy- at (ao,b ). Consequently, the functions are 



continuous there, or 

(8) lim (f>k(a + Aa, b + Ab) = cj>k{a ,b ) , A; = 1,2 

The limit of a complex function is assumed to be a well-understood notion. If we rewrite 
(Q) in the form of complex pairs using po — (ao,£>o) and Ap = (Aa, Ab) we get an 
equality, whose validity follows from (Js[) . In other words, tp is continuous at po. 

In Chapter 4 we shall also need the following formula. Let (f> : (a,b) — > <fi(a,b) be a 
complex function that is holomorphic in a domain GCB. The continuity of ^ and ^ 
means that at any point (a, b) € G we can write 

(9) 0(a + Aa, b + Ab) - <f>(a, b) = ^ ■ Aa + ^ ■ A6 + e(Aa, A6) 

where e(Aa, A6) is a complex-valued function that tends to more rapidly than Aa and 
Ab tend to 0. 

3.2 The exponential function 

With the help of the complex exponential, sine and cosine functions we define the bi- 
complex exponential function by 

(10) e p = (e a -cos6,e a -sino) ,p=(a,b) 

The component functions </>i(a, b) — e a ■ cosb and </>2(a, 6) = e a • sin b are obviously 
holomorphic in both a and b. Definition ( |To| ) is of the same form as the corresponding 
definition of the complex exponential function rendered as a pair of reals: 

(11) e a — (e x ■ cosy, e x ■ sin y) ,a—(x,y) 

With b = (2, u) the complex functions cos 6 and sin& have the real pair representations: 

cos& = (cos z coshu, — sinz sinhit) 
sin& = (sinz cosh u, cos zsinhu) 

Together with (|ll]) they enable us to rewrite (|l0|) as a quadruple in R 4 . For p = (a, b) = 

(2, y, z, u) we get 

(12) e p = ^1,^2,^3,^4) 

= ( e x (cos y cos z cosh u + sin y sin z sinh u) , 
e x (— cos y sin z sinh it + sin y cos z cosh it) , 
e x (cos y sin z cosh u — sin y cos z sinh it) , 
e 31 (cos y cos z sinh u + sin y sin z cosh u) ) 
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Turning our attention to the properties of e p we first observe that for p = (a, 0) we get 
e ( a :0) _ (e<^ o). The bicomplex function thus generalizes the right complex function. 

Definition (|l(]) also gives us the periodicity of e p . In the complex plane e a has the 
period im2n, to = ±1,±2, . . . , while cos b and sin b have the period n2n, n = ±1,±2, . . .. 
The period w of e p is therefore (im, n)2n, formally: 

(13) e p+w =e p torw = (im,n)2vr 

to = 0,±1,±2,... 
n = 0,±1,±2,... 

Setting p = yields as special case 

(14) e (im,n)2ir _ x ,m = 0,±l,±2,...andn = 0,±l,±2,... 

The familiar addition formula e a+b = e a ■ e b satisfied by the complex exponential function 
is valid in bicomplex space, too, provided complex multiplication is replaced by bicomplex 
multiplication: 

(15) e p+r = e p Oe r , p = (a, b) and r = (c, d) 

For r = —p this formula becomes e p e~ p — 1, which shows that e p never gets the value 
zero. In Section 3.6 we shall obtain the stronger result that e p never gets a singular 
value. 

3.3 The hyperbolic and trigonometric functions 

With the exponential function at our disposal we are now in a position to define the 
hyperbolic functions cosh and sinh in B in the same way as it is done in C: 

(16) coshp = ^(e p + e- p ) ,p=(a,b) 

(17) sinhp = l( e p -e- p ) 

The use of (^) in these formulas yields after straightforward calculation the complex 
pair representations 

(18) coshp — (cosh a ■ cos b, sinh a ■ sin b) ,p~(a,b) 

(19) sinhp = (sinh a ■ cos b, cosh a ■ sin b) 

Note that for p — (a, 0) the functions reduce to their complex analogues. Moreover, due 
to their definitions they have the same period as e p . 
Bicomplex division (2.^9|) gives us the tanh-function 

sinh v , , , 

(20 tanhp = , p = (a, b) 

coshp 
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By applying @, (|T|) and (2.g) in this formula we get after some manipulation 

tanh a tan 6 



(21) tanhp = 



(cos b) 2 + (tanh a ■ sin 6) 2 ' (cosh a) 2 + (sinh a • tan b)' 



To define the bicomplex cosine and sine functions we take the corresponding complex 
functions as starting point: 

(22) cosa = ^(e ia + e- ia ) , a = x + iy 

(23) sina = ^-(e w - e~ m ) 

Employing real pairs to express complex numbers we have a = (x, y), i = (0, 1) and 

ia = (0, 1) • (x,y) = {~y,x) 

Transferring this to B we replace a = (x, y) by p — (a, b) and i = (0, 1) by j = (0, 1) to 
obtain 

(24) j Qp — (0,1) ©(a, b) = (-6, a) 

We choose the expression j p as the bicomplex correspondence of the complex ia and 
generalize ( p^ ) and (^3|) by 

(25) cosp = i( e J©P + e - J0p ) ,P=M) 

(26) sinp = (e J0p - e~ 3&p ) 
Their C 2 -representations are 

(27) cosp = (cos a • cosh b, — sin a • sinh &) ,p=(a,b) 

(28) sinp = (sin a • cosh b, cos a • sinh b) 

Because cosa and sina both have the period m27r, m = ±1, ±2, . . . , while coshfr and 
sinhfe have the period in2ir, n = ±1,±2, . . ., the period of cosp and sinp is (m,m)27r. 
The B- and (^-representations of the bicomplex tan-function are: 

sin p 

(29) tann = , p = (a, b) 

cosp 

(30) t ( tana tanh b 

\ (cosh 6) 2 + (tana ■ sinh 6) 2 ' (cosa) 2 + (sina • tanh b) 2 , 

In the definitions of the preceding functions we observe in every case that they are the 
same as the definitions of their complex counterparts but with complex notions replaced 
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by bicomplex ones. This isomorphism manifests itself in the properties of the functions, 
too. For example, the familiar addition formulas of cosh, sinh, cos and sin can directly 
be taken over to bicomplex space: 

cosh(p + r) — coshp cosh r + sinhp sinh r 

sinh(p + r) = sinhp cosh r + coshp sinh r 

cos(p + r) = cosp0 cosr — sinp0 sinr 

sin(p + r) = sin p cos v -\- cos p sin v 

We also note that the functions are interconnected through the formulas: 

cosp = cosh(jQp) 

sinp — — sinh(j © p) 

i 

3.4 Polynomials 

Bicomplex polynomials are formed as sums of bicomplex integer powers p n multiplied by 
bicomplex constants. 

As the simplest polynomial we regard the constant function 

(31) — po , po a bicomplex constant 
The bicomplex identity function is 

(32) I(p) = (a,b) ,p = M) 
and the bicomplex power function 

(33) F(p) = (a, b) n , p — (a, b) and n an integer > 



Here the exponentiation should be performed according to (2.27). We have in particular 

p° = (l,0) 

p 2 = {a 2 -b 2 ,2ab) 

The components of the complex pairs of these equations are examples of so-called har- 
monic polynomials of a and b [Q. These polynomials are also recursively definable: if we 
denote them by G n (a,b) and H n (a,b) we have 

(34) p n = (G„{a,b),H n (a,b)) , p = (a, b) and n > 
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From the identities p° — (1, 0) and p n+1 = pQ p n we then derive the recursive scheme: 

G (a,b) = 1 , H o (a,b) = 

G„+i(a, b) = a ■ G n (a, b) - b ■ H n (a, b) 

H n+ i(a,b) = a- H n (a,b) + b ■ G n (a,b) 

The general form of an n th degree bicomplex polynomial with argument p — (a, b) is thus 

(35) Q(p) =d n @p n + dn-i p" _1 + . . . + di p + d Q 

where dk stands for a bicomplex constant. 

Similarly, one can define bicomplex power series of the form (Sfc : k > : dk &p k )- 
Their meaningfulness depends as usual on their convergence. 

3.5 The quotient function 

Let the bicomplex functions ip and 9 be given by 

^(P) = (<t>i(a,b),<f>2(a,b)) , p = (a,b) 
6(p) = {uj l {a,b),uj 2 {a,b)) 
For the inversion Q^ 1 we first have on application of (2.|38|) 
<:-!(„ < Wl 



The function is well-defined in the domain of B where u)\ + ui\ ^ 0, that is where 
CN((uji, L02)) 7^ 0. A bicomplex inversion can have rather more zeros in the denominator 
than the corresponding complex function, a point to remember when one attempts to 
generalize e.g. Cauchy's integral formula. 

The inversion of the identity function (^) is a special case of (||(]) : 

The quotient function ip/8 = ip # _1 has the form 

Jp_ _ f 4>i ■ u>i + <h. ■ w 2 4>2 ■ u)i - 4>i ■ u 2 
^ > a ~ \ , ,2 1 , ,2 ' , ,2 1 , ,2 

An important category of quotient functions consists of rational functions of the type 
(39) R ^ = W) 

where Ql(p) and Q2(p) are bicomplex polynomials without common factors. Borrowing 
terminology from complex function theory we say that the zeros of Q2(p) are the poles 
of R{p). 
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3.6 The logarithm function 

We shall define the bicomplex logarithm function as the inverse of the bicomplex expo- 
nential function. The definition will be based on the following result. 

Theorem 3.1 The bicomplex number q = (7 + i5, e + irj) can equivalently be rendered in 
the form 

(40) q = (v, 0) (cos w, sinw) 

where v and w are complex numbers, if and only if q is nonsingular or equal to zero. 

Proof. It is without loss of generality possible to express v and w in terms of the real 
numbers r, y, z, u so that 

(41) v = re ly , w = z + iu 

Our task is therefore to investigate under what conditions the equation (7 + i5,e + in) 
= (re ly ,0) (cos(z + m),sin(z + iu)) has a solution in the unknowns r 7 y,z,u. The 
right-hand side is equal to 

-re iy ■ (e i( - z+i "' ) + e -' i ( z + iu )^ —re iy ■ (e^ z+iu ^ — e _ ^ z+itl ^)^ 

and by rewriting this expression using the addition formula as well as other basic prop- 
erties of the complex exponential function, one can show that the equation is equivalent 
to the following system of equations: 

^r(e~ u cos(y + z) + e u cos(y — z)) = 7 
^r(e~ u sin(y + z) + e u sin(y - z)) = 5 
ir(e~ u sin(y + z) — e u sin(y — z)) = e 

-r(—e~ u cos(y + z) + e u cos(y — z)) = n 

If q = 0, i.e. 7 = (5 = e = ?7 = 0, these equations are solved by r = 0, in which case 
the values of y, z and u are indeterminate. If q 7^ and singular, i.e. of the form 
(7 + i5, — S + i-y) or (7 + iS, S — i-f), the equations admit no solution (since they imply 
either sin(y + z) — A cos(y + z) — or sin(y -z) = 0A cos(y — z) = 0). However, if q 
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r = 


(K-L) 1 ' 4 


K 


= (i + v) 2 + (s- 


£ ) 2 




y = 


\{M + N) 


L 


= (i-v) 2 + ($ + 






z = 


\{N-M) 


M 




- 717T 


, n = V n = 1 


u = 


Mnl^l 

4 1 L 1 


N 


= arctan (|±|) H 


- 717T 


, n = V n = 1 



is nonsingular they are solved by: 



(42) 



The nonsingularity of g ensures that K and -L are both nonzero, thus making u well- 
defined. When M is evaluated the principal value in the interval (—§,§) should be 
chosen for the arctan-function. The integer n should be set to or 1 depending on 
whether the denominator 7 + 77 is positive or negative. If the denominator is zero M = ^ 
or M — depending on whether the numerator 6 — e is positive or negative. The same 
applies, mutatis mutandis, to the evaluation of N. 

The theorem follows from the existence of the solution (fl2"|). 

□ 

In the bicomplex number representation q = (v,0) (cost/;, sin u;) we call v the scale 
factor and w the complex argument of q. We shall use the notations: 

v — s f(<l) {scale factor of q} 

w — ca(q) {complex argument of q} 

The function ca(q) is multi- valued, because its real part is not uniquely determined; for a 
given q its infinitely many values differ from each other by multiples of 2ir. It is possible 
to identify a principal branch of ca(q), denoted by Ca(q), that is characterized by the 
fact that its real part is in a certain interval, [0, 27r) say. Formally: 

w = Ca(q) = w — ca(q) A < Re(w) < 2n 

The complex exponential function gets every nonzero value in C, hence if q ^ in formula 
( |40| ) v can be replaced by e a for some a & C. This gives q = (e a - cosu>, e a ■ sin w ) or the 
bicomplex exponential function. With Theorem 3.2 we then deduce that this function 
gets every nonsingular value and no singular value. 

We denote the bicomplex logarithm function by blog and define it as the solution of 
the equation e p = q: 

blog(g) = p = e p = q for p, q G B 

Clearly, blog(g) is only defined for nonsingular q. If we set p = (a, b) and lean on (|40| ) 
and the definition of the exponential function, e p = q becomes 

(e°,0) (cos&, sin 6) = (v, 0) (cos w, sin w) 
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This equation is satisfied if and only if 

a = log(w) 

b = w + ?i27r , n an integer 

where log stands for the complex logarithm function. By making use of the scale factor 
and complex argument of q — see ( ff3| ) — we further obtain: 

a = log(sf(q)) , q nonsingular 

b = ca(q) + n2ir , n an integer 

Thus, we have derived the following expression for blog: 

(44) blog(q) = (log(s/(g)) , ca(q) + n27r) , q nonsingular and n an integer 

Observe the isomorphism between this function and its complex counterpart log(i>) = 
In \v\ + i arg(w). 

The blog-function is multi- valued, because both the log- and the ca-function are multi- 
valued. If p is a bicomplex logarithm of the nonsingular bicomplex number q, then so are 
p+ (im, n)2ir for all integers m and n. The principal branch of blog, which we denote by 
Blog, is a function, however, and consists by definition of the principal branches of the 
component functions: 

(45) Blog(g) = (Log(s/(g)) , Ca(q)) , q nonsingular 

The principal branch Log of the complex logarithm function is defined in the normal 
manner. 

The bicomplex logarithm function allows us to extend the definition of the power 
function q r to bicomplex exponents: 

(46) q r — e r0Blog ^) for q, r G B and q nonsingular 

Example. Just for fun we shall evaluate (0, The bicomplex number (0, i) is 

nonsingular, hence with the help of (|4l]) and (^2|) we may render it in the form (|40|): 

(e l 2, 0)0 (cos-, sin-) 

This yields 

Blog((0,i)) = (Log(e^),|) = (z|,|) 

and further 

(0, i)( Q '') = e (O,i)0Blog((O,i)) _ e (-if ,-f) _ ,q^\ 

□ 
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We are also able to define the inverse functions of the trigonometric functions in the same 
way as it is done in C. For example, the inverse cosine is obtained as the solution of 

cosp= I( e J©P + e ^ P) =q for p,g£ B 

Solving for e J0p and taking the bicomplex logarithm yields the bicomplcx version of a 
familiar formula: 

arccos q = —j © blog(g ± \J q 2 — 1) 
It is of course required that blog's argument q ± ^/ q 2 — 1 is nonsingular. 
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4 Differentiation of bicomplex functions 



We give the derivative of a bicomplex function a definition that is isomorphic to the 
corresponding definition in C and because, in addition, the bicomplex operations are 
isomorphic to the complex ones, we may expect that the normal rules for computing the 
derivative of polynomials, sums, products, quotients, etc. remain valid in B. However, 
in defining the bicomplex derivative we have to reckon with the fact that the bicomplex 
numbers do not form a field in the ordinary sense due to the lack of inverses of the 
singular numbers. The possibility to restrict the discussion to nonsingular numbers only 
must be excluded, since they do not form a subfield in B: if a and b are two nonsingular 
bicomplex numbers, then a + b need not in general be nonsingular, although a © b is so 
(see formula (|35|)). In spite of this, the problems related to singularities can be solved. 

We shall obtain B-, C 2 - and R 4 -representations for the bicomplex derivative. The 
C 2 -representation is particularly useful and will give us the bicomplex version of the 
Cauchy-Riemann equations, which every differentiable, bicomplex function has to satisfy. 
We shall show that such a function possesses derivatives of all orders. 

4.1 Definition of the derivative 

Following the standard procedure one is tempted to define the derivative of the bicomplex 
function ip at a point p as the limit 

V U ' Ap^O Ap 

A problem then arises, because the fraction (ip(p + Ap) — ip(p))/Ap is not defined for 
each singular Ap in the neighbourhood where the limit is taken. An obvious remedy is 
to confine limit taking only to nonsingular Ap, which demands that we apply a weaker 
formula than (3.||) for the limit concept, namely 

(1) (VAp, e : CN(Ap) / 0Af>0: (38 : S > : ||Ap|| < S => 

ip(p + Ap) - ip(p) 



Ap 



-V'Cp) 



<e)) 



The dummies satisfy Ap £ B and e, 6 G R. Recall from Section 2.5 that Ap is nonsingular 
if CN(Ap) 0. The above limit condition is denoted by 

lim feMijW 

Ap^O Ap 

Formula (Q) is an instance of the following limit condition that pertains to bicomplex 
quotients in general: 
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(2) (Vp, e : CN(0(p)) ^ A e > : (3<5 : c5 >0 : 

\\p-po\KS- m 

The corresponding limes notation is 



0(p) 



w 



<e)) 



w = km — - 
p->-po 0(p) 

CN(e(p))=£0 

After these preliminaries we can state our first version of the definition of the bicomplex 
derivative: 

Definition 4.1 Let ip be a bicomplex function whose domain of definition contains a 
neighbourhood of the point p. The derivative of tp at the point p is defined by the 
equation 

(3) iP>(p)= lim tk+M^tM 

Ap^O Ap 
CW(Ap)^0 

provided this limit exists. □ 

A bicomplex function that has a derivative at the point p is said to be differentiable or 
holomorphic at p. If the function is holomorphic at all points of a domain G C B it is 
said to be holomorphic in G. The derivative ip'ip) is as usual also written 

Remark. By a domain we understand in the sequel a connected, open set of points in 
the bicomplex space (occasionally in the complex plane). □ 

Definition 4.1 could at least formally be strengthened by the requirement that for every 
singular Ap such that p + Ap is in the neighbourhood of p in question, the fraction 
Q(Ag) = (ip(p+ Aq) — ip(p))/Aq should approach a finite limit w as Aq approaches Ap, 
formally: 

w= lim tk+M^M ,CN(Ap) = Q 

Aq^Ap Aq 

GN(Aq)=iO 

Note that Aq must be nonsingular. Note also that different values of w can be associated 
with different values of Ap. The existence of such a limit means that it is possible to 
remove the singularity of f2 at each Ap by defining Q(Ap) — w, which makes the condition 
CN(Ap) 7^ superfluous in (||). Nevertheless, this strengthening just postpones the 
moment when one has to compute a limit in a neighbourhood consisting of nonsingular 
points only. 

Example. Let ip be the bicomplex square function ip(p) = p 2 , which yields 

ijjjp + Ap) - ip(p) _ 2p Q Ap + (Ap) 2 
Ap~ ~ ~Ap 
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The limit value 2p + Ap of the right-hand side exists for all singular Ap. Of these Ap = 
gives the derivative. □ 



Yet another alternative is to introduce the derivative of "0 by a definition of Frechet- 
type: if there exists a bicomplex number K and a bicomplex function e that satisfy the 
equations limAp^o e (Ap) = 0, e(0) = and 

(4) ip(p + Ap) - i(j(p) = KQ Ap + e(Ap)0 Ap 

then the derivative ip'(p) is equal to K. Thanks to the absence of denominators in 
(^), the problems related to singular bicomplex numbers are at first sight circumvented. 
However, the identification of a suitable function e presupposes in general the existence 
of removable singularities of the fraction (ip(p + Ap) — ip(p))/Ap. The equivalence of this 
definition and Definition 4.1 has been demonstrated by Price [f23f , who refers to (||) as 
the strong Stolz condition. 

Independently of which of the above definitions is applied, one obtains by employing 
standard techniques the normal rules for computing the derivatives of sums, products 
and quotients of functions: 

(V + 0)'(p) = ^'( P ) + 0'(p) 

(V©0)'(p) = iP'(p)®6(p)+Tp(p)Q6'(p) 

|)'w - m^&zmim .cm m) ^ 

Moreover, for the bicomplex composite function (tjjo6)(p) = tp{6{p)) the chain rule holds: 

(V°0)'(p)=V>'(0(p))o0'(p) 

The derivative of the power function is given by the normal formula 

,r\ ^ n-1 

(5) ~r = n P 

dp 

provided n is a positive integer and p any bicomplex number or n is a negative integer and 
p a nonsingular bicomplex number. In the next section this derivative will be generalized. 

4.2 Complex pair formulas for the derivative 

We shall seek formulas for the derivative of a bicomplex function tp having the complex 
pair representation 

(6) ^(p) = (cj>i(a,b),(t>2(a,b)) ,p = (a,b) 
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The calculations are essentially the same as the ones that lead to formulas (3.||) and 
(3.^) for the complex derivative, but due to the possibility of singular denominators it is 
necessary to go through them in detail. 

We investigate what limits definition (||) yields if we let Ap approach in two ways 
so that Ap = (Aa, 0) or Ap = (0, Ab). Then, a key observation is that, provided Aa ^ 
and Ab ^ 0, Ap is never singular. In other words, we do not need to worry about 
singularities in the subsequent formulas. 

On account of (2.^) the inverse of Ap is 

(7) J- = (^-,0) ifAp=(Aa,0) 



Ap V Aa 

(8) ± = (0,-£) if Ap = (0, Ab) 

With respect to the case Ap = (Aa, 0) we get using (||) and (Q) 

V>(p + Ap) — ij)(p) _ { fa(a + Aa,b) — fa(a,b) fa(a + Aa, b) - fa{a, b) 



Ap V Aa Aa 

As Ap and Aa approach simultaneously, the components of the right-hand side ap- 

^ and ^ 

aa da 



proach the complex derivatives and The bicomplex derivative therefore takes 



the form 

(9) -dp-^y-da-^) >P=^ b ) 

If we deal with the second case Ap = (0, Ab) in the same way we get with (||) and (||) 

4){p + Ap) - 4)(p) _ ( fa{a,b + Ab)- fa{a,b) fa (a, b + Ab) - fa (a, b) s 



Ap \ Ab ' Ab 

Letting Ap and Ab approach results in 
, 1fV . dip (dfa dfa\ 



dp \ db db 

This formula must be equivalent to (^]), a requirement whose consequences will be exam- 
ined shortly. 

We are now equipped to find more derivatives of the elementary functions. Applied 
to the exponential function e p — (e a ■ cos b, e a ■ sin b) both (|) and @ give 

de p 

(11) ^- = e p 

dp 

In order to find the derivative of the principal branch of the bicomplex logarithm function 
(3.[l5|) we assume that the four bicomplex numbers p, Ap, q, Aq satisfy q — e p A q+Aq = 
e p+Ap , thereby implying that q and q + Aq are nonsingular. This yields 

Blog(g + Ag) - Blog(g) Ap 

Ag ~ e p+Ap _ e p 
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As Ap tends to 0, the right-hand side tends to l/e p . At the same time Aq tends to 
and the left-hand side becomes the derivative of Blog: 

fw\ dB\og(g) 1 . 

(12) ; = - , q nonsingular 

dq q 



With the help of (jllj) , (J12j) and the chain rule we obtain the derivative of the generalized 
power function (3.[46|) 

dq r 

(13) — — = r q r , q nonsingular 

dq 

For integer r this formula is equivalent to (^|). One should keep in mind that if r is a 
positive integer, the requirement that q is nonsingular can be dropped, a fact of particular 
interest in the development of power series, for instance. 

We also mention that the complex pair representations (3jl^), (3 jlSj), (3.27), (3.p8[) of 



the hyperbolic and trigonometric cosine and sine functions yield the expected derivatives: 

dcoshp 



dp 

dsinhp 
dp 

dcosp 
dp 

dsinp 



sinhp 
coshp 
— smp 
cosp 



dp 

4.3 The bicomplex differentiability condition 

The requirement that the two complex pair formulas ([)]) and (|Io|), which we found for 
the bicomplex derivative, should be equivalent is met if and only if the equations 

(14) ^1 = ^ 
{ ' da db 

(15) -da~ = —db 

hold in the domain G C B where ip(p) = {<pi(a, b), (f>2{a, b)) is differentiable. These 
equations represent the bicomplex differentiability condition and will be referred to as 
the bicomplex Cauchy-Riemann equations (abbr. bicomplex CR- equations). They are the 
complexification of the CR-equations (3.||). 

The validity of ( |l4| ) and (|l5| ) presupposes that the partial derivatives and 
k = 1, 2, exist, i.e. that the functions are holomorphic in a and b. That this is indeed 
only a necessary condition for -0 to be holomorphic in p is illustrated by the following 
example. 
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Example. Let 6 be given by 

0(p) = (a 2 ,b 2 ) ,p=(a,b) 

The functions <pi(a,b) = a 2 and 02 (a, 6) = b 2 are holomorphic in a and b, but 9 is not 
holomorphic in p, because 4>\ and 4>2 do not fulfill the bicomplex CR-equations. These 
are instead fulfilled by the components of the square function 

if>(p) =p 2 = (a 2 -b 2 ,2ab) 

which thus is holomorphic in p. □ 

Note that all bicomplex functions introduced in Chapter 3 fulfill the bicomplex CR- 
equations. 

From the above we conclude: 

Lemma 4.1 If the bicomplex function ip(p) = b), ^2(1, b)), p— (a, b), is holomor- 

phic in the domain GCB, then the complex functions (f>\ and <p2 are holomorphic in G 
and satisfy the bicomplex CR-equations (JJ^[j-(|J^) there. □ 

We shall demonstrate that the converse also holds: 

Lemma 4.2 // the complex functions </>i(a,b) and </> 2 (a, b) are holomorphic in the do- 
main GCB and satisfy the bicomplex CR-equations JT^J-dJ^) there, then the bicomplex 
function ^p{j>) = (0i(a, b), </> 2 (a, b)), p — (a,b), is holomorphic in G. 

Proof. The holomorphism of k — 1,2, means that they have partial derivatives 
and that are continuous in G. At an arbitrary point p = (a, b) of G we may therefore 
by virtue of (3|^) write 

<^(a + Aa,b + Ab)-0i(a,b) = ^ ■ Aa + ^ • Ab + e 1 

oa ob 

2 (a + Aa,b + Ab)-0 2 (a,b) = ^ • Aa + ^ • Ab + e 2 

oa ob 

where e\ and e 2 are complex-valued functions of Aa and Ab such that they tend to 
more rapidly than (Aa, Ab) in the sense that (ei(Aa, Ab), 0) / (Aa, Ab) — > and 
(0, e 2 (Aa, Ab)) / (Aa, Ab) -> as (Aa, Ab) -> 0. If we set Ap = {Aa, Ab) and make use 
of 

tp(p + Ap) - i(i(p) = ((f>i(a + Aa,b+ Ab) - cj)i(a,b) , (t>2(a + Aa,b+ Ab) - (j} 2 (a,b)) 
together with (|lj) and ( |l5[ ) we are now able to derive 

^( P + Ap)-^(p) = >f)0Ap + ( £l , £2 ) 
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Dividing both sides by Ap and splitting the last term into two terms gives 
0(p + Ap) - if,(p) fdfa 90 2 \ , (ei,0) , (0,e 2 ) 



Ap \ da ' da J Ap Ap 

The assumptions about ei and €2 imply that the last two terms tend to as Ap tends 
to 0, which shows that 4>'(p) exists and equals 

9^1 902 
da ' da 

Here, the computation of limits is in each case restricted to nonsingular Ap in accordance 

with (0) and (§). 

□ 

Combining the preceding lemmas results in: 

Theorem 4.1 Let ip(p) = (0i(a, 6), 02(a, 6)), p = (a, b), be a bicomplex function. Then 
ip is holomorphic in the domain G C B if and only if the complex functions 0i(a, b) and 
2 (a, 6) are holomorphic in G as well as fulfill the bicomplex CR-equations ^L^j-^lh]) 
there. □ 

The bicomplex CR-equations are essential to the theory under consideration. They can 
be combined into a single equation, namely 

(16, .'-<*>> 

where 2* and S£ denote 

dj)_ _ / 90i 902 \ 9^ _ / 90i 902 

[ ' da~ \ da ' 9a / ' 96 ~ I 96 ' 96 



Formula (16) is the generalization of || + i • ^ = satisfied by a complex holomorphic 



function. 

4.4 Higher order derivatives 

Next we investigate the existence of higher order derivatives of ip. Application of (j^) 
times to ip gives 

(18) ^=U^'9^J 

Alternatively, the n th derivative can be computed by repeated application of (|Tc|): 
M £ = ((-l)L.j^,(-l)L*j^) ,forodd^l 



(20) ^ = f ( -i)Lfj^i, M) m^ 

V 7 d P n y ' 96" ' V ; 96" J 



, for even n > 
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For the proof of the subsequent theorem we note that ( |l8| ) can be rewritten 

•">-<> 

The left-hand side is conveniently abbreviated by 

(22) ^ = p , n > 

It is understood that tp^ = ip. 

Theorem 4.2 A bicomplex function ip{p) = (<fii(a, b), $2(0-, &)), V = ( a ?^)> * 5 
morphic in a domain G C B possesses derivatives of all orders in G, each one holomor- 
phic in G. 

Proof. Let the n th derivative of ip be given by ([l8]) and designated by ip( n \ According to 
Lemma 4.2 is holomorphic in G if and ^ a t 2 are holomorphic in G and satisfy 
the bicomplex CR-equations. 

By Lemma 4.1 the holomorphism of -0 in G implies that 0i and 02 are holomorphic in 
G. From complex analysis we know that a complex holomorphic function has derivatives 
of all orders, which by themselves are holomorphic. It follows that the holomorphism of 
0i and 2 is inherited by ^ and 

To prove that 0( n ) satisfies the bicomplex CR-equations we apply ([l6]), instead of 
and establish 

using mathematical induction. For n = the formula is identical with (|l6|), which holds 
on account of 0's holomorphism. For the validation of the step we assume ( p3|) and 
calculate: 

g^(n+i) _ ^ «9V> ( ™ +1) _ 

{ ^(n+i) = fig^, ^ @ and 



da \ da J db \ da 

{induction hypothesis} 

d ( dip( n) \ d (dtb^\ 

^{-J Q ^) + l G d-b{—) =0 

{taking derivatives in different order, see below} 

d fdtp^\ d fdip^\ 

- 3& db[^a-) +3(3 ¥ b [^a-) =0 
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= {identity} 
true 

In the penultimate step we changed the order of the operators and ^ when applied 
to ip( n \ The interchange is allowed because consists of holomorphic component 
functions and 

d_ fduA. _ d_ /dw\ 
da\db J db\da 



holds for an arbitrary complex function lu : (a, b) — ► u>(a, b) that is holomorphic in both 

a and b. 

□ 

4.5 R 4 -representations of the bicomplex derivative 

In order to get another handle to bicomplex derivatives we rewrite (^|) and ( |l0| ) in R 4 . 
Assume that the full representation of the function tp is 

ip(p) = {<t>i{a,b),(j)2{a,b)) 

(<p 1 (a,b),(j)2(a,b)) = (ipi(x, y, z, u) + iip2(x, y, z, u) , 
ip 3 (x, y, z, u) + iip 4 (x, y, z, u)) 

p = (a, 6) , a = x + iy and b — z + iu 

On account of (3.||) and (3.[|) the partial derivatives and are given by: 

d x ipi + id x ip 2 



da 
da 

db 
db 



dytp2 - idylpl 

d z ipi + id z ip 2 
dui>2 - idui>i 



Analogous formulas hold for and Substitution of the complex derivatives into 
(||) and ( ^p| ) yields four R 4 -representations of 

(24) ^ = (d x ipi + id x ip 2 , d x ip 3 + id x tpA) 
dp 

(25) ^ = {d y ij)2 - idyip! , dyip 4 - id y ip 3 ) 
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(26) ^j- = {d z ip3 + id g ip4 , -d x ipi - id z ip 2 ) 
dp 

(27) ^j?- = (9„-04 - id u ip 3 , -d u ^% + id u ^x) 
dp 

Formulas @ and (§§) were obtained from (g), @ and @ from @. 

Example. The R 4 -representation of the square function ip(p) = p 2 is 

i/)(p) = (x 2 — y 2 — z 2 + u 2 + i(2xy — 2zu) , 2.tz — 2yu + i(2yz + 2ccu)) 
Each of (|2^)-(p7|) yields the same derivative 

^- = (2.t + i2j/, 2z + i2u) 
dp 

This means that = 2p, as expected. □ 



The formulas fl24|)-([27|) should be equivalent. The equivalence of (|2J) and (|25|), in par- 
ticular, implies the simultaneous validity of the equations 





= dy1p2 




= -d y 1pl 




= dyljj4 


d x ip 4 





They express that ipi + iip2 and tp3 + wpi are holomorphic with respect to a = x + iy. By 
the same token, the equivalence of d2g) and fl27|) amounts to the simultaneous validity of 



5*^3 = 9„V4 

d z i>4 = -d u ip 3 
—d z ipi = -<3 U ^>2 

which express that ipi + iip 2 and ip3 + ^4 are holomorphic with respect to b — z + iu. 
The additional equivalence of ( pi[ ) and ( |26| ) means that the equations 

d x ip3 = -d z ipi 
dxi>4 = -d z ip 2 
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all hold and validate the bicomplex differentiability condition (fi^)-(|l5|). Thus, in a 
straightforward way formulas (p4j)-(p7|) exhibit both the holomorphism of i/j and its 
component functions ifii + i?/>2 and ^3 + iipi- 
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5 Bicomplex integration 



In this chapter we define the concept of a bicomplex integral as a line integral. We inquire 
under what conditions it is path independent using an argument analogous to the one 
given by Ahlfors for complex functions. By integration we show that a holomorphic 
bicomplex function can be written as a Taylor series. To prove Cauchy's theorem in B 
we lean on a generalization of Green's theorem. To express Cauchy's integral formula in 
B we first derive the bicomplex twining number, which generalizes the complex winding 
number. We are then able to prove the formula with the help of a Taylor expansion. 

5.1 Definition of the line integral 

We consider the integration of the bicomplex function 

(1) = (Ma,b),Ma,b)) 

((j> 1 (a,b),<p 2 (a,b)) = (^lO^y, z,u) +iip 2 {x,y, z,u), 
ip 3 (x, y, z, u) + iip 4 (x, y, z, u)) 

p = (a,b) , a = x + iy and b = z + iu 

For the time being ^>'s domain of definition is left anonymous. We assume that (j>\ and 
<p2 are holomorphic in a and b, thereby ensuring that ip is continuous as demonstrated 
in Section 3.1. 

The basic bicomplex integral is essentially isomorphic to the complex integral. It is to 
be understood as a line integral that is evaluated with respect to some four-dimensional 
curve T in B. More specifically, the concept to be defined is 



J V'(p) © dp , dp — (da, db) 



Henceforth, we shall choose the curve T so that it is piecewise continuously differcntiable 
in B and has the parametric equation 

(2) r : p = p{t) , p(t) = (a(t),b(t)) for r < t < s 

The argument t is real. One can look at T as a curve made up of two component curves 
7i and -f 2 in C 

(3) r = ( 7l , 72 ) 

whose parametric equations are 

71 : a = a(t) , a(t) — x(t) + iy(t) for r < t < s 

72 : b = b(t) , b(t) = z(t) + iu(t) for r < t < s 
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As definition of the line integral of ip(p) extended over the curve T we then take 



iP(p)Qdp= / ip(p(i))®p'(t)dt 



(4) 



Because ip is continuous, ip(p(t)) at the right-hand side is also continuous. If p'(t) is 
discontinuous at some points the integration has to be performed in subintervals of [r, s] 
in the normal manner. 

The left-hand side of (@) can be reformulated in the following two ways. Firstly, 
bicomplex multiplication of tf> = (<f>\, (f>%) and dp = (da, db) yields 



(5) J tjj(p)Qdp — I J 4>i ■ da — 02 ■ db, J 4> 2 ■ da + 4>\ ■ db 

Secondly, since the differential dp is expressible as 

dp= (da,0)+jQ(db,0) ,j=(0,l) 



we have 
(6) 



ip(p) © dp 



tp(p) © (da, 0) + j ip(p) © (db, 0) 



The bicomplex differentials (da, 0) and (db, 0) have the same properties as the complex 
differentials da and db. Formulas (|J) and (||) will be needed in the next two sections. 

We shall also need integrals of the type J r ip(p) \\dp\\ involving the norm (2.f42f). With 
T specified by (0) we define 



m\\d P \\ = f Wt))ib'(*)ii 

J r 



dt 



The special case ip(p) = 1 gives us a suitable definition of the length of V: 



(8) 



\dp\\ = length of T 



By regarding the integral (|J) as the limit of a Riemann sum we obtain by applying the 
triangle inequality 



(9) 



ip(p) dp 



< 



\Hp)\\ ■ \\M\ 



If the integral at the left-hand side has a finite value there exists a real constant M such 
that 11^(^)11 ^ M f° r sull p on T. We therefore get the estimate 



(10) 



tp(p) dp 



< M-L 



where L stands for the length of T. 

Formulas (|7|)-(|l0|) are direct generalizations of the corresponding ones in complex 
analysis p8[. 
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5.2 Path independence and the primitive function 

We ask ourselves under what condition the bicomplex line integral is independent of the 
path of integration in a domain GCB. With respect to definition (^) this means that if 
r is free to vary in G the integral depends only on the end points of T. In particular, if 
T is a closed curve the integral reduces to zero. 

Our reasoning will be based on formula The fact that its right-hand side is of 
the form 

J uji (da, 0) + lu 2 Q (db, 0) 
gives relevance to the following theorem. 

Theorem 5.1 Let u>i and 0J2 be two bicomplex junctions of p = (a, b) such that they are 
continuous in the domain GCB. Then the line integral J r ^i (da, 0) + u>2 (db, 0), 
defined in G, depends only on the end points of the curve T if and only if there exists a 
function fl : (a, 6) — > 0(a, b) such that ^ = u>\ and ^ = uj2- 

Proof. Omitted, because it is entirely isomorphic to the proof of the corresponding 

theorem for complex functions as given by Ahlfors pp. 106-107. 

□ 

We focus on the evaluation of the right-hand side of (](]) in a domain GCB. According 
to the theorem above the integral depends only on the end points of T if and only if there 
exists a function "J : p — > ^(p), where p = (a, b), such that 

(id 

(12) «=ie#) 



Multiplying (|12|) by j and adding it to ([Ll]) yields 

9* a* 

oa ob 

But this is the bicomplex differentiability condition (4|l^), which means that '5 is a 
holomorphic function. Formula ([ll]) further tells us that ip is the derivative of ']/. (Note 
that ^ — §7 on account of (4.fjl|), and that ip, too, is holomorphic on account of 
Theorem 4.2.) Hence, we have established: 

Theorem 5.2 Let ip be a bicomplex function of the form ([!]) that is continuous in a 
domain GCB. Then for an arbitrary curve T in G the integral § T ip(p) dp depends 
only on the end points of T if and only if ip is the derivative of a holomorphic function 
* in G. □ 
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The function '5 is uniquely determined up to an additive constant and may be chosen as 
primitive function (integral function) of ip. It enables one to evaluate a line integral of 
tp extended over a curve with end points p\ and P2 according to 

fPI 



i>{p)Odp = [*(p)]£ = *(p2)-*(pi) 



Thanks to path independence the omission of any reference to the actual curve is justified 
here. 

5.3 Taylor series 

We shall show that a holomorphic bicomplex function can be developed into a Taylor 
series. Our starting point is the formula for computing the derivative of the product of 
two bicomplex functions: 

{0Q<p)'(p) = e'(p)Q<p(p) + 6(p)Q<p'(p) 

In a domain G where both 9 and <p are holomorphic the function (OQip)' certainly has a 
primitive function, or 9 ip. As a result, if we subject both sides to the integral operator 
f£ o and transfer a term between them we obtain the bicomplex version of the formula 
for integration by parts: 

r 9{p) q<p'(p) Qd P = \6{ P ) <p( P )}p o - f V(p) © ipip) © dp 

jpa 

It is possible to compute the integrals along any curve in G starting at po and ending at 
P- 

Let the function tp be holomorphic in G so that it possesses derivatives of all orders 
there due to Theorem 4.2. Then, because integration by parts is available, we can expand 
tp into the following Taylor polynomial and remainder term about p a : 

(13) V(P) = (Vfc:0<fc<n-1: (p - Po ) k ) +R n (p,Po) 



■p 



Ru(p,po) - t ~~ttt / Tp (n) (q) (p - q) n ~ L © dq 

(n - !) ! Jpo 

In order to estimate the remainder term R n (p, po) we study its behaviour in a po-centered 
ball that lies in G and is specified by \\q — po\\ < \\p — po\\- We assume that the integral is 
evaluated along a curve whose length \\dq\\ equals L. The derivative ip^ n \q) satisfies 
II (q)\\ < M ior & suitable M. Furthermore, using the triangle inequality, in particular, 
we deduce: 

lb-g|| = \\(p-Po) - (q ~Po)\\ < lb - Po|| + Ik -Poll < 2\\p-po\\ 
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Hence, we estimate 

„ , \ i M ■ L ■ (2\\p- Poll)™ -1 

(14) \\Rn(p,P0)\\ < 

Because ||i?„(p,po)| approaches as n — > oo the infinite Taylor series converges. 
5.4 Cauchy's theorem in bicomplex space 

For the purpose of generalizing Cauchy's theorem for a bicomplex function we recall one 
of its formulations in the complex plane. 

Theorem 5.3 Let A be a simply connected domain in C. If the complex function f(a) = 
g(x, y) + ih(x, y), a — x + iy, is holomorphic in A, then 



f(a) -da = Q 



for any closed curve 7 in A. □ 



The domain A being simply connected means that it is without "holes" or, more precisely, 
that its complement with respect to the extended plane is connected. It is well-known 
that the theorem is provable with the help of Green's theorem, provided one assumes that 
g{x, y) and h(x, y) possess continuous partial derivatives in x and y §, §|. Although the 
proof is based on unnecessarily strong assumptions, we take it as starting point because 
conducting it in B clarifies important concepts. In R 2 Green's theorem is written 

(15) I g dx + hdy = — — I dx dy 



J 1 Js( 7 ) \dx dy/ 

The left-hand side is a line integral extended over a closed curve 7, the right-hand side a 
surface integral extended over $(7), which is the domain of the x-y plane bounded by 7. 

Transferring the above to B first demands that we get some basic understanding of 
complex line and surface integrals of the form 

(16) / cp-da + uj-db 



(17) / <j)-da-db 

where 

• r is a closed curve in B with the parametric equation (^). 

• S(T) is an arbitrary surface in B bounded by T. 

• (j> and to are complex functions of a and 6, i.e. of the type B — > C (or C 2 — ► C, if 



one so wants). 
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Integrals like ( |16[ ) appeared already in formula (|5|). For the present purpose it suffices 
to define the integrals ([l6]) and ( ]l7| ) by reducing them to real line and surface integrals, 
which we consider well-understood. Let <f> and u> more specihcally be given by 

(j)(a,b) = ^(x,y,z,u) + i £ 2 (x,y, z,u) 

uj(a,b) = r]i{x, y,z,u) + i r] 2 (x, y,z,u) 

where £1, £2, Vi, V2 are continuous functions of the type R 4 — > R. With da = dx + idy 
and db — dz + idu we obtain 

(18) J <p ■ da + ui ■ db = J (£1 dx — £2 dy + 771 dz — r/2 du) 

+ i J (£2 dx + £,idy + r/ 2 dz + rji du) 



(19) / (f) ■ da - db = / (£1 dxdz — £1 dydu — £2 dxdu — £2 dydz) 

JS{T) JS(T) 

+ i (£1 dxdu + £1 dyefo: + £2 cteiiz — £2 dydu) 
Js(r) 

We generalize Green's theorem ( |T5| ) by substituting complex functions and variables for 
real ones as well as the four-dimensional curve T for the two-dimensional 7. We also have 
to require that the functions are holomorphic. 

Theorem 5.4 Let (f>\ and 4>2 be two complex functions of a and b such that they are 
holomorphic in a domain G C B. Also let S(T) be an orientable surface in G bounded 
by the closed curve T . Then 

(20) / ct> 1 -da + ct> 2 -db= [ f^l-^l) -da-db 
Jr Js(T) \9a db J 

Proof. Uses the generalized Stokes theorem, see Appendix. 
□ 

With this result the proof of Cauchy's theorem in B comes out in few lines. 

Theorem 5.5 Let the bicomplex function ip(p) = ((f>i(a, b), (f>2(a, b)), p— (a,b), be holo- 
morphic in a domain GCB. Then 



ip(p) Qdp = 

r 

for any closed curve T that is the boundary of an orientable surface S(T) in G. 
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Proof. 



ip(p) © dp 

{(!)} 



i • da — (f>2 ■ db, / (f>2 ■ da + <f>i ■ db 
i Jr 



{S(T) is an orientable surface in G; ( ppj ) twice} 

-^a--^b)- da - d ^j s{T) \-da-~^b)- da - dh ) 
{ip fulfills the bicomplex CR-equations (4jl|)-(4.|l|) in G} 







□ 



Example. By way of illustration we shall evaluate the integral of e p along the closed 
curve 

T: p(t) = (Re lt , Re lt ) , < t < 2ir 

where R is a positive real constant. On the x-y- and z-u-planes T is projected as two 
origo-centered circles of radius R. The component curves of T have the parametric 
equations 

71 : a(t) = Re lt , < t < 2ir 

72 : b(t) = Re lt 

The outcome of the following calculation complies with Cauchy's theorem. 
e p dp 

{bicomplex line integral, < t < 2w} 
e p(t) &p'(t)dt 
{definition of T and e p } 



r 



2tt 

(e Re%1 • cos(i?e 1 *) , e Re%% • sin(i?e 1 *)) (iRe lt , iRe lt ) dt 



{multiplication, distribution properties} 
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2tt 

(e Relt ■ cos(Re lt ) - e Re " ■ sin{Re lt )) ■ iRe lt dt , 

(e Reit • sin(i?e lt ) + e ReU • cos(i?e lt )) • iRe lt dtj 
{definition of the curves 71 and 72} 
(^J (e aW - cos a(t) - e a ^ • sin a(t)) • a'(t) dt , 

jf (e 6 <*>-sin&(i) + e h ^ ■ cos b{t)) ■ b'(t) dt^j 
{definition of line integral in C} 
I / (e a - cos a — e a - sin a) • da , / (e b - sinfe + e b ■ cos 6) • dfe 
{71 and 72 are closed curves; Cauchy's theorem in C} 





□ 



5.5 The twining number 

In the complex plane the so-called winding number u(-y, do) of a closed curve 7 with 
respect to a point ao is defined by 



*(7, ao) = f da , ao $ 7 

2m J a — a 



We call the corresponding bicomplex concept the twining number. The twining number 
of a closed curve T with respect to a point p in B is denoted by v(T,p ) and defined by 

(21) v(T,p Q ) — f — — — , p — po nonsingular for all p e T 

2ttj JrP-Po 

The integral has a finite value only if p — po is nonsingular at every point p on T. The 
value cannot always be found by application of Cauchy's theorem, and may therefore 
differ from zero. The application of the theorem demands that T is the boundary of 
some orientable surface S(r) so that p — p n is nonsingular for all p € S(T), but such a 
surface does not exist, if T can be regarded as enclosing p . The notion of curve enclosure 
is more difficult to grasp intuitively in the bicomplex space than in the complex plane, 
and the twining number is in fact the right device for making it precise: 

Definition 5.1 The curve T is said to enclose the point p in the bicomplex space if 

w(r,po) ^0. □ 
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Example. We shall compute the twining number v(T,po) when V is given by 

r : p(t) =p + {e u ■ cost, e i{ • sin t) , < t < 2tt 

Notice as an aid to the geometric interpretation of this curve that all its points are at a 
distance of 1 from po- 

We first observe that p(t) —po = (e lt ■ cos t, e lt ■ sin t) — e^**'*' . Because the exponential 
function gets only nonsingular values, p(t) — po is nonsingular for all t. We also have: 

1 / cost sint N 



This gives: 



p(t) — po \ e lt e l 
p'{t) = (i ■ e lt ■ cost - e lt ■ sint, i ■ e u ■ sint + e lt ■ cost) 



dp 



rP-Po 

{line integral, < t < 2ir} 

&p'{t) dt 



p(t)~Po 

{formulas above and bicomplex algebra} 

(i, 1) dt 



= {integration} 
(i27r,2vr) 

Thus, the twining number ( pT| ) becomes v(T,po) = (1, — i). It differs from 0, which means 
that F does not bound any surface S(T) that would have the property that l/{p — po) is 
holomorphic at all points of S(T). For instance, take the surface 

5(r) : p(h,t) =p + {h-e lt - cost, h ■ e lt ■ smt) 

where the real parameters h and t are in the intervals [0, 1] and [0, 27r], respectively. It is 
bounded by T, but contains the point p(0, t) = po, where l/(p — po) is not holomorphic. 
□ 

In order to find a general expression for the twining number we evaluate the right-hand 
side of (^l]) assuming that T's parametric equation is (||) and that p — p is nonsingular 
at all points p £ T. We introduce the function 

(22) F{t) =fP^)_ dt 

Jr PW ~ PO 



43 



F(s) is the value of the integral in (pl|). F is continuous on the interval [r, s] and except 
at those points where p'(t) is discontinuous it has the derivative 

**(*)- P ' {t) 



p(t) - Po 

As a result the derivative of the expression e~ F ^ (p(t) — po) vanishes except at a finite 
number of points. The expression itself therefore reduces to a constant, namely p(r) — po- 
Accordingly, we obtain 

e F(t) _ Pit) - Po 
p(r) - po 

For t = s we have e F W = 1, since V is closed and p(s) = p(r). To fulfill this identity 
F(s) must due to (3.[l4|) be the period of the exponential function. Consequently, setting 
t = s in (p3) yields first 

/ ^ - - — dt — (im,n)2ir , m and n are integers 

Jr P(t) - Po 

and then, on omission of the parametrization 

/ — — — = (im,n)2TT , m and n are integers 

Jr P - Po 

Hence, (|lj) becomes 

(23) v(T,po) = (to, —in) , m and n are integers 
For n = the twining number reduces to the complex winding number. 

5.6 Cauchy's integral formula in bicomplex space 

We are now ready to tackle the bicomplex generalization of Cauchy's integral formula. 

Theorem 5.6 Let tp be a function, p Q a point, T a closed curve and G a domain in the 
bicomplex space and assume that they satisfy the following conditions: 

a) ip is holomorphic in G. 

b) r is contained in G and the difference p — po is nonsingular for all p G Y . 

c) There exists an orientable surface S(T) contained in G and bounded by the 
curve r so that po & S(T) A v(T,po) — holds andp — po is nonsingular for 
all p £ SiT), or po <G >5(r) A v{T,po) ^ holds andp — po is nonsingular for 
allpeS(T)-{po}. 

Then, we have 

(24) ^(Po) © v(T,p ) = -L / dp 

2ttj Jr P ~ Po 
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Proof. Our argumentation is essentially the same as the one given by Nevanlinna and 
Paatero for the complex case pl| , pp. 133-135. It is primarily based on conditions a) 
and c), condition b) being a necessary condition of c). 

Condition c) gives us two options with regard to the surface S(T). If it does not pass 
through p 0l the function 4>{p)/{p — Po) is holomorphic at all p <E S(T). Formula ( ]24| ) 
then holds, since v(T,po) = 0, and the integral at the right-hand side is equal to due 
to Cauchy's theorem. 

If S(T) passes through p and v(T,p ) ^ 0, i.e. T encloses poi we consider the function 

ip(p) - ip(p ) 



F(p) = 



if P ^ Po 

p-po 

ip'ipo) tfp = Po 

By virtue of conditions a) and c) F is holomorphic at all points of the punctured surface 
S(T) — {po}, while the situation at p$ remains open. Our aim is to apply Cauchy's 
theorem to F and for this purpose we must show that F is holomorphic at po, too. 

The function ip is holomorphic at j»0i hence it has a Taylor expansion of the form ( |l3| ) 
about this point: 

V>(P) = ^Po)+^'(Po)©(p-Po) + ^"(Po)©(p-Po) 2 + 

«V> '"(Po) © (p - Po) 3 + Ri{p,Po) 
6 



We get 



f(p) = tM^M = ^(p ) + V(p )o(p-p o) + 

P -po 2 



1 , III, n „ / VI , Ri(.P,Po) 

-xV (Po)0(p-Po) H 

6 p- po 



and further, because F(pq) — ip'(po): 



F(p) - F(p ) 1 1 Ra{p,Po) 

= oV (Po) + ftp (Po) (p - Po) + t 

p-Po 2 6 (p-poT 

The inequality ( |l4| ) implies that i? 4 (p,p )/(P ~ Po) 2 tends to as p — > po- We conclude 
that F'(p ) = ^tp"(po), or that F is holomorphic at po- 

Since F is holomorphic on the whole surface S(T) it is legitimate to use Cauchy's 
theorem 

F(p) Qdp = 



which yields 



o = / tM_zlM Qdp = f ^PL Qdp - Hpo)@ f d P 



P-Po JrP-Po JrP-Po 
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Formula (24) follows if we apply (21) and substitute v(T,po) 2ivj for the last integral. 
□ 



Assuming that v(T,po) ^ in ( pif ) and that conditions a) - c) of the preceding theorem 
otherwise also hold, we may select a particular curve T so that v(T,po) = (0, — i) in 
agreement with (p3[) to obtain 



(25) 4( po ) = *Q[MPLQdp 

2m Jr P-Po 

We conclude that for such a curve the common form of Cauchy's integral formula is valid 
in B, too. 

Example. Given the closed curve 

r : p{t) = {e i \e lt ) , < t < 2ir 

we shall compute the value 

1 f e p 
W = — © / — dp 
2m J r p 

By making the identifications ip(p) = e p , po — and G — B we see that the conditions 
a)-c) of Theorem 5.6 are met. Especially concerning c) we note that v(T, 0) = (0, — i) 
and that all points p S V are nonsingular (since e %t ^ ie lt A e lt ^ — ie %t for all t £ [0, 2ir]). 
Moreover, the surface 

S(T) : p(h, t) = (h-e lt , h- e lt ) , < t < 2n and < h < 1 

passes through the point p(0,t) = and all points p € S(T) — {0} are nonsingular. 
Therefore, formula (^5|) is applicable and yields W = e° = 1. We arrive at this result by 
direct calculation, too: 

p 

-Qdp 
r P 

{bicomplex line integral, < t < 2ir} 

^ e p(t) 

, — —Qp'(t)dt 

o P(t) 

{definitions of e p and T; bicomplex inverse} 
C\e elt ■ cos(e Ji ) , e e " • sin(e lt )) (j^ , © ( ? e li , ze Jt ) d£ 

{multiplication, distribution of the integral} 



' - cos ( eit ) t At M /- 2 "e elt -sin(^) „ 



2tt 

ie rt dt, / — ■-■ic"rfl 
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= {71 : a(t) = e lt and 72 : b(t) = e lt are the components of T} 

/ n -a'(t)dt , / ,, N w -b / (f)dt 

Uo o(*) W Jo &(*) / 

= {definition of line integral in C} 

/ f e a -cosa , /" e b -sin& \ 
/ da ' / — j; d& 

= {Cauchy's integral formula in C} 

(27ri,0) 

The result is W = 1, as it should. 
□ 



47 



6 Bicomplex harmonic analysis 



Our last topic is the relationship between the bicomplex functions and the Laplace equa- 
tion. It is well-known that the complex holomorphic functions satisfy the two-dimensional 
Laplace equation, hence one may expect that the bicomplex holomorphic functions sat- 
isfy its four-dimensional version. We shall show that this is indeed true for eight classes 
of bicomplex holomorphic functions. One of them consists of the functions that fulfill 
the bicomplex Cauchy-Riemann equations (4jl^)-(4.|l5|), the others of similar functions 
that are partly or entirely formed with the conjugates of complex holomorphic functions. 
These functions quite naturally emerge from the factorization of the Laplace equation 
in its various forms. We start by factorizing the two-dimensional equation, since it will 
provide us with important concepts for the more general cases. 

6.1 The Cauchy-Riemann operator 

With / a complex function depending on the real variables x and y, the two-dimensional 
Laplace equation is written 

d 2 f d 2 f n 
— - H = 

dx 2 dy 2 

or, using abbreviated differential operators 

(1) d 2 J + d 2 y f = 

The two-dimensional Laplace operator A2 = d 2 + d 2 admits the factorization 

(2) A 2 = {d X + idy) ■ (0 X - idy) 

which for the complex number a = x + iy gives us the occasion to introduce the so-called 
Cauchy-Riemann operator D a Q : 

(3) D a — d x + idy , a = x + iy 

Note that D* = d x — id y . 

We now let the complex function / in (1) more specifically be given by 

(4) /(a) = g{x,y) +ih(x,y) ,a = x + iy 

where g and h are functions of the type R 2 — > R such that they have continuous partial 
derivatives in x and y. From (2) and (3) we obtain the relation 

(5) A 2 f(a) = <= (D a ■ f(a) = 0) V (D* ■ f(a) = 0) 

by means of which the Laplace equation of the consequent can be solved by solving either 
disjunct of the antecedent. We first evaluate D a ■ f(a) = 0: 



48 



D a ■ f(a) = 

= {definition of D a and /, omitting arguments} 

{d x + id y ) ■ (g + ih)=0 
= {complex multiplication} 

d x g - d y h + i(d x h + d y g) = 
= {complex algebra} 

{d x g - dyh = 0) A (d x h + d y g = 0) 

The first disjunct of the antecedent of (||) obviously resolves into the Cauchy-Riemann 
equations 

(6) d x g = dyh , d x h = -d y g 

An analogous treatment of the second disjunct D* ■ f(a) = shows that it is equivalent 
to 

(7) d x g = -dyh , d x h = d y g 

We call these the conjugate Cauchy-Riemann equations. The name is justified by the 
fact that they are solved by the conjugates of the holomorphic complex functions. These 
are henceforth referred to as the conjugate holomorphic functions. 

Example. The exponential function g + ih = e x cosy + ie x suiy is holomorphic and 
satisfies @. Its conjugate (g + ih)* = e x cosy + i(—e x sin y) satisfies (0). □ 

From the above we deduce the following relationship between the Cauchy-Riemann op- 
erator and the holomorphic/conjugate holomorphic complex functions: 

(8) D a ■ f(a) = = / is a holomorphic function 

(9) D* ■ f(a) = = / is a conjugate holomorphic function 

The application of the Cauchy-Riemann operator or its conjugate to a complex function 
thus reveals whether the function is differentiable or not. 

6.2 The derivatives of a complex function 

It is a noteworthy fact that given a complex number d with non-zero real and imaginary 
parts any complex number a can be expressed in the form a = ad + f3d* , where a and (3 
are real numbers. Because d and its conjugate d* span the entire complex plane in this 
way, one could argue that in the theory of complex functions complex-valued variables 
and functions ought to be treated on an equal footing with their conjugates. This should 
especially apply to the derivatives of the functions. With the complex function / specified 
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by (||) it is possible to distinguish between four derivatives: 

df_ dt_ df_ df*_ 

da ' da* ' da* ' da 

We first recall the definition of the first of these derivatives at the point uq = xq + iyo 

(10) ± = lim /(a) = /(ao) 

da \a.=a a— >a a — Of) 

Letting a approach ao along the real and imaginary axis separately leads to the formulas 
(the evaluation at the point ao now omitted) 

(11) £ = d x g + id x h 



(12) -L = d y h-id y 



4f 

da 

which are equivalent if the CR-equations (||) hold. 

Next we consider the derivative of the conjugate function /* with respect to the 
conjugate variable a* . We define 

d£_ = Um r(a)-f*(a ) 

da* |a*=aS a*-*a* a* — Uq 

But if the limit at the right-hand side of ( |l0| ) exists we have 

um /*w-/*M = ( lim /(<*) - fMY 



a*-^-a* a* — Oq ya— >a a — &o 

from which it follows that 

(13) ^* - ( ^ 

da* \da 

Conjugation of ([ll]) and ( |l2|) consequently results in two formulas for 4^: 

df* 

(14) — = d x g - id x h 

df* 

(15) — = dyh + id y g 

The equivalence of these equations is also implied by the CR-equations. 

Then we turn our attention to the derivative of / with respect to a* . Its definition is 

£ = Um f (Q) ~ /(ao) 

da \a*=a* a*^a* a* — <2q 

The right-hand side can be rewritten 

,. 9{x, y) - g(x 0) yo) + i(h(x, y) - h(x 0) y )) 
lim ; — 

x - iy ^ x-iy-(x - iyo) 

xo—iyo 
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Letting a* approach a*, along the real and imaginary axis amounts to setting y — yo and 
x = x , respectively, in this formula. Taking the two limits separately yields 

(16) — = d x g + id x h 

(17) ^ = -d y h + id y g 

These equations are equivalent if the conjugate CR-equations (]?]) hold, i.e. if / is conju- 
gate holomorphic. 

For the derivative finally, we deduce with the same technique that lead to ( [l3| ) 

(18) dt = (d[_ 
V ' da \da* 



Thus, by conjugating ( |16| ) and (|T7|) we get 

df* 

(19) — = d x g~id x h 

da 

df* 

(20) -4- = -d v h-id v g 

da y y 

whose equivalence likewise follows from the conjugate CR-equations. 
Summarizing these results, we have: 

Formulas @/@ and @/@ are pairwise equivalent if / is a holomorph 



complex function 
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Formulas (16)/fll7j) and (tt9j)/(E0j) are pairwise equivalent if / is a conjugate holo- 



morphic complex function. 

The preceding argument admittedly contains some redundancy, because we have ob- 
tained four pairs of equivalent formulas, although we have only two types of functions, 
holomorphic and conjugate holomorphic ones. The reason is that we have designated a 
conjugate holomorphic function by both / and /*. The duplicity could of course easily 
be removed, but experience has shown that one needs both notations and therefore all 
the formulas, too. 

It is clear that in the complex plane the conjugate holomorphic functions must possess 
properties that are very similar to those of the holomorphic functions. The conjugate 
CR-equations (0) play the same role for the former functions as the CR-equations play 
for the latter; in practice this means that if a relation R(x, y) has been shown to hold for 
an holomorphic function / of the form (|j) the relation R(x, —y) will hold for /*. Due to 
this kind of close similarity it can be regarded as an effort devoid of interest to actually 
develop a theory of conjugate holomorphic functions. What we shall mainly need from it 
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in the sequel is that a conjugate holomorphic function / of the form has a continuous 
derivative computable with ( |l6| ) or (]l7|). Alternatively, if the function is designated 
by /* we can apply or ([[5]) to find the derivative. 



6.3 Factorization of the four-dimensional Laplace equation 

We return to the bicomplex space and consider once more a function tp of the type 

(21) = (4>i(a,b),h(a,b)) 

((j>i(a,b),cf)2(a,b)) = (ipi(x, y, z, u) + itp 2 (x, y, z, u), 
ip 3 (x, y, z, u) + itpi(x, y, z, u)) 

p = (a,b) , a = x + iy and b = z + iu 

We want ip to fulfill the four-dimensional Laplace equation 

(22) Efcl> + d 2 v ^ + dli> + 9> = 
which we shall solve by factorization. Using the Laplace-operator 

we render (S3) in the form 

(23) A 4 x^ = 

where x denotes quaternionic multiplication (2.^). The A^operator is to be regarded as 
a scalar quaternion, hence the multiplication in ( ]23] ) merely applies it to each component 
ip k of ij> — see (2.©. 

We now make the crucial observation that the Laplace operator admits the factoriza- 
tion 

(24) A 4 = (d x + id y ,d z + id u ) x (d x - id y , -d z - id u ) 
If we introduce the new operator L q by 

(25) L q = (d x + id y ,d z + id u ) 
the factorization becomes 

(26) A 4 = L q x L* 

The Cauchy-Riemann operators D a = d x + id y and Db = d z + id u give us 

(27) L q = (D a ,D b ) , L* = (£)*,- D 6 ) 
and 
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(28) A 4 - (D a ,D b ) x (D* a ,-D b ) 

With formulas (|24|), (Jji^) and ( ^8| ) we have factorized the four-dimensional Laplace op- 
erator in R 4 , B, and C 2 , respectively. 

Because L q x L* = L* x L q we deduce from ( p3| ) and ( pG] ) 

(29) A 4 x?/> = (L 9 x %j) = 0) V (£* x ip = 0) 

by means of which we can solve the Laplace equation by solving either disjunct of the 
antecedent: 

(30) L q x tp = 

(31) L*xip = 

We focus on the first of these. With ( p7| ) we first rewrite it in C 2 : 

L q X tp = 
= {L q = (D a ,D b ),ip = (4> 1 ,4> 2 )} 

(D a ,D b ) x #a)=0 
= {quaternionic multiplication} 

(A, • 0! - A, • 0* , D t • # + D a ■ fc) = 
= {equality of quaternions} 

(D a ■ 4>i - D b ■ q)* 2 = 0) A (D b ■ <j}{ + D a ■ <j) 2 = 0) 



Thus, (|3C|) is equivalent to the simultaneous validity of the equations 

(32) ZV0 a - D b -cj>* 2 

(33) D b ■ (f>* = -D a ■ 4> 2 

Substituting D a = d x + id y , D b = d z + id u and 4>\ = ipi + iip? , 4>2 — "03 + iipi m these 
we obtain with complex algebra the corresponding R 4 -representation: 

(34) d x il>i - dyip 2 - d z ip 3 - duipi = 

(35) d x ip2 + dyipi + d z ip4 - d u ip 3 = 

(36) d x ip 3 - d y ip4 + d z ipi + d u ip 2 = 

(37) d x ip 4 + dyip 3 - d z ip 2 + d u ip! = 

This partial differential equation was derived in 1935 by Fueter ju^. As mentioned in 
the Introduction his starting point was the formula 

(38) d x -0 + i x d y ip + j x d z ip + k x d u ijj = 
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which we get by substituting the vector representation of L qi or L q = d x + d y i + d z j + d u k, 
in L q x tp = 0. 

Equations (|30|), (j32|)-([33"|), (|3^)-(p7|), and ( |38|) are four equivalent ways of expressing 
the same condition. We shall refer to them as the Fueter equations. 

The treatment of equation ([H]) , L* x ip — 0, which also enables one to solve the Laplace 
equation, is analogous with the first case. Written in the form (D*, —D b ) x (<^x, <fe) = 
— see (^lj) and ( p7| ) — it becomes equivalent to the simultaneous validity of 

(39) D* a -4>i = -D b -4>1 

(40) D b -4>\ = D* a -fc 

We call these the conjugate Fueter equations. Expanded in R 4 they look like: 

(41) 3^1+^2 + ^3 + ^4 = 

(42) d x il> 2 - d y tpi - d z ip 4 + d u ip 3 = 

(43) d x il> 3 + d y ip4 - d z ipi - d u il> 2 = 

(44) 3^4-^3 + ^2-^1 = 

This partial differential equation was derived in 1929 by Lanczos Jig ]. The fourth, equiv- 
alent formulation of the conjugate Fueter equations is found by substitution of the vector- 
form L* — d x — d y i — d z j — d u k in L* x ip — 0, which becomes 

(45) d x tp — i x d y if) — j x d z ip — k x d u tp — 
6.4 Regular and conjugate regular functions 

In the solutions of the Fueter equations we identify a new type of bicomplex functions, 
hence we define: 

Definition 6.1 A bicomplex function is said to be regular in a domain of the bicomplex 
space if it satisfies the Fueter equations in this domain. □ 

To further characterize the regular functions we prove: 

Theorem 6.1 Let the bicomplex function if) be defined in the domain G C B so that 
ip(p) = (0i (a, 6), 4>2(a, b)) for p = (a, b). If the functions <p\ and 4>2 are holomorphic in a 
and conjugate holomorphic in b at all points (a, b) G G, then i\) is regular in G. 

Proof. The assumptions about </>i and <j)2 imply on account of (||) and (JoJ) that 

D a -<p k -0 , fc= 1,2 

D* b -<t> k = , fc = l,2 
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hold in G. Conjugating the second formula yields -D& • <j)* k = 0, which combined with the 

first validates the Fueter equations (|32|) and (33). Hence, ip is regular in G. 

□ 



An example of a regular function is 

(46) E(p) = (e Q -cos*6,e a -sin*6) ,p = (o,6) 

The complex pair components e°- cos* b and e°- sin* 6 are functions of the variables a and 
b such that they are holomorphic in a and conjugate holomorphic in b, as required. 

Two-argument, complex functions that are holomorphic in the one argument and 
conjugate holomorphic in the other are evidently of special interest and deserve an own 
term: 

Definition 6.2 A two-argument, complex function </> : (a, b) — > <fi(a,b), defined in a 
domain G C B, is said to be coanalytic in G if it is holomorphic in a and conjugate 
holomorphic in b at all points (o, b) £ G. □ 

Due to Theorem 6.1 the coanalyticity of the components of a bicomplex function ip is a 
sufficient condition for ^'s regularity. 

It is a fundamental property of any conjugate holomorphic complex function /* that 
f*{b) = f(b*)- The right-hand side of (||^) can therefore be rewritten (e a -cos6*, e a -sin&*), 
which we also obtain by replacing the argument p — (a, b) of the exponential function 
(3.0) by q= (a, b*): 

(47) e q = (e a -cos&*,e Q -sin6*) , q = (a,b*) 

With the same change of argument we derive the regular counterparts of the other ele- 
mentary bicomplex functions of Chapter 3, e.g. the regular identity function and power 
function — see (3. [32]) and (3.f33|): 

(48) I(q) = (a,b*) ,q={a,b*) 

(49) F(q) = (a, b*) n , q = (a, b*) and n an integer > 

and the regular hyperbolic and trigonometric functions — see (3.^8|)-(3.p^) and (3]2 

(50) coshq = (cosh a • cos 6*, sinha • sin b*) ,q=(a,b*) 

(51) sinhq = (sinha • cos b* , cosh a • sin b*) 

(52) cosq = (cos a • cosh 6*, — sin a • sinho*) 

(53) sinq = (sin a • cosh&*,cosa • sinh&*) 
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Note that all the component functions of the complex pairs at the right-hand sides are 
coanalytic in B. 

The holomorphism of one-argument complex functions is preserved by complex ad- 
dition, multiplication and division as well as by functional composition. The conjugate 
holomorphism of one-argument functions is similarly preserved. Being a combination of 
these properties, the coanalyticity of two-argument complex functions is thus preserved 
by the aforementioned operations. 

Example. The functions a ■ cos b* and b* ■ sin a are coanalytic in C 2 , hence a ■ cos b* + 
b* ■ sin a and (a ■ cos&*) • (b* ■ sin a) are coanalytic in C 2 . □ 

We are now able to prove that the regular functions form an algebra under the bicomplex 
operations and that, in addition, functional composition preserves the regularity property. 

Theorem 6.2 Let the bicomplex functions ip : q —> ip(q) and 9 : q — > 9{q), where q = 
(a, b*), be regular in a domain GCB, and assume that they can be represented by complex 
pairs of coanalytic functions. Then their bicomplex sum ip + 9, product ip®9 and quotient 
tp/0 are regular in G, the last only at those points qo where CN(9{qo)) 0. Furthermore, 
the composite function (ip o 9)(q) = ip(9(q)) is regular in G, provided the domain of ip is 
in the range of 9. 

Proof. We verify the regularity of the product ip 9, in particular. According to the 
assumption ip and 9 have representations of the form 

^(l) = (<t>t(a,b*),<fa{a,b*)) ,q=(a,b*) 

9{q) = (uj 1 (a,b*),aj 2 (a,b*)) 

where the functions ef>i,<f>2, <+>i and u>2 are coanalytic in G. Bicomplex multiplication gives 

1p 9 = (01 • LOi - 02 ' ^2 , 4>2 ■ + 4>l ■ W2) 

Because complex addition and multiplication preserve coanalyticity, the component func- 
tions 0i • ojx — 4>2 ' ^2 and 02 • k>i + 01 " ^2 are coanalytic in G, thereby making ip 9 
regular. The regularity of ip + 9, tp/9 and ip o 9 is verified by analogous arguments. 

□ 

Next we focus on the conjugate Fueter equations, especially in their form (^)-(^o|). A 
bicomplex function satisfying these equations is said to be conjugate regular. We prove 
a sufficient condition for a function to be conjugate regular (compare Theorem 6.1). 

Theorem 6.3 Let the bicomplex function ip be defined in the domain G C B so that 
ip(p) = (0i(a, 6), 02(a, 6)) for p = (a, b). If the functions 0i and <p 2 are conjugate holo- 
morphic in both a and b at all points (a, b) G G then ip is conjugate regular in G. 
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Proof. The assumption about <j>\ and <f>2 implies on account of (0) that 

D* a -ct> k = , A; =1,2 

I? b *-0 fc =O , A; =1,2 

hold in G. The first formula and the conjugation of the second then together validate 
□ 

If we take the quaternionic conjugate of the regular identity function ( |48| ) we get I*(q) = 
(a* , —6*) or, by a change of argument 

(54) I(q*) = (a*,-b*) ,q* = (a*,-b*) 

We regard this as the conjugate regular identity function. In the same way, by applying 
the other elementary bicomplex functions of Chapter 3 to the argument q* = (a*, —b*) 
we derive their conjugate regular counterparts, for instance: 



(55) 


cosh q* 


= (cosh a* 


• cos b* . 


— sinh a* • sin b*) 


,q* = (a*,-b*) 


(56) 


sinh q* 


= (sinh a* 


■ cos b*, 


— cosh a* • sin £>*) 




(57) 


cos q* 


= (cos a* • 


cosh b*. 


sin a* • sinh b*) 




(58) 


sin q* 


= (sin a* • 


cosh b*, 


— cos a* • sinh b*) 





The component functions of the complex pairs at the right-hand sides are conjugate 
holomorphic in both a and b. 

Analogously to the regular functions the conjugate regular ones form an algebra under 
bicomplex operations. 

6.5 New classes of differentiable functions 

By considering the derivatives of the regular and conjugate regular functions and of simi- 
lar functions we shall obtain altogether eight classes of holomorphic bicomplex functions. 

Several authors have studied the derivative of regular functions in the context of 
quaternionic function theory (see e.g. f28|). A general observation has been that it is 
unproductive to define the derivative of a regular quaternionic function ip at the point p 
by applying definition (2.^l|) for the quaternionic inverse and taking the limit 

lim [ty(p + Ap) - il>(p)) x ((Ap) T -1)] 

Ap— >0 

because for p = x+yi + zj +uk and Ap = Ax+Ayi + Azj + Auk the class of differentiable 
functions will become too restricted to be of interest. (The class will contain only first- 
degree polynomials. Nothing is gained by using left-division instead of right-division in 
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the above expression.) The remedy is to employ bicomplex division and to insist on 
taking the limit of a fraction that is regular in Fueter's sense. The latter is achieved by 
replacing p and Ap by the regular 

q = (a, b*) = (x + iy, z — iu) 

Aq = (Aa, Ab*) = (Ax + iAy, Az - iAu) 

Hence, the derivative of a holomorphic regular function is introduced as follows. 

Definition 6.3 Let ip be a regular bicomplex function whose domain of definition con- 
tains a neighbourhood of the point q = (a, 6*). The derivative of tp at the point q is 
defined by the equation 

A g ^o Aq 

CN(Aq)^0 

provided this limit exists, when it is computed for nonsingular Aq = (Aa, Ab*). □ 

Apart from the function argument this definition is the same as Definition 4.1 for the 
derivative of an "ordinary" bicomplex function of the variable p = (a, b). This means that 
if one makes the appropriate change of argument, the expected formulas for computing 
the derivative of a sum, product and quotient of bicomplex functions as well as the 
chain rule are retained for the regular holomorphic functions. Furthermore, if tp has the 
complex pair representation tp(q) = ((px(oL, b*), <fa((i, &*)) its derivative is also given by 
the following equivalent formulas, analogous to (4.^) and (4.[lO|) : 

( 59 ) ~7Z = I ' "n7 ) , q={a,b) 



dq \ da da 

1 J dq \db* ' db* 

Application of either formula to the regular exponential function e q = (e a -cos b* , e a -sin b*), 
for example, yields the same result or 

de q 

— = e q 
dq 

The derivatives of the other elementary regular functions are likewise found with the 
normal formulas. 

Remark. In the application of (pOf) it is helpful to remember ( p^|) 

da* \ da 
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satisfied by a conjugate holomorphic complex function /*. It enables one to compute the 
derivative of e.g. cos b* according to: 



dcosb* d(cosb)* ( dcosb 



db* db* V db 

□ 



(— sin b)* = — sin b* 



The equivalence of (59) and (p0[) presupposes the simultaneous validity of the equations 



(61) -fr = W 

1 ' da db* 

They express the differentiability condition of the regular functions and are termed the 
regular Cauchy-Riemann equations. 

Example. For q = (a,b*) the function 

#z) = (a 2 ,(&*) 2 ) 

is regular but not holomorphic, because it does not fulfill the above equations. In contrast, 
the square function 

iP{q)=q 2 = {a 2 -{b*Y,2ab*) 
is both regular and holomorphic. □ 

It is instructive to rewrite (^9|) and (|60|) in R 4 assuming that <p% = -01 + iip2 an d <p2 — 



ip3 + iip4- Using formulas (|l l|) — (pL 2[) and (pL6|)— ( 17 ) of Section 6.2 we get two formulas 
for each of the derivatives 3^., §g£ and §|j. Substituted into (§|) and @ they 
furnish four R 4 -representations of 

(63) ^j- = {d x tpi + id x ip 2 , d x ip 3 + id x ipi) 
dq 

(64) ^ = {dyip 2 - idyipi , dyipi - id y ip 3 ) 

(65) ^j- = {d z ip 3 +id z ip4, -d z ipx - id z ip 2 ) 
dq 

(66) ^j- = {-d u ip4 + id u ip 3 , d u ip 2 - id u ipi) 
dq 

The first two of these were obtained from ([39]), the last two from (|60|). 
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Example. For q = (a, b*) = (x + iy, z — iu) the R 4 -representation of the regular square 
function ip(q) = q 2 is 

tp(q) = (x 2 -y 2 - z 2 +u 2 + i(2xy + 2zu) , 2xz + 2yu + i(2yz - 2xu)) 

Each of (|6^)-(|66]) yields the same derivative or 

-i- = (2x + i2y,2z-i2u) 
dq 

This means that ^ = 2q, as expected. □ 

The central concepts pertaining to the differentiation of the conjugate regular functions 
are analogous. It was shown in the previous section that the application of bicomplex 
functions to the argument q* — (a*,—b*) yields the corresponding conjugate regular 
functions, as illustrated by (|5^)-(|58|). The same argument is used in the definition of the 
derivative of such a function. 

Let tp(q*) = (0i (a*, —b*), 4>2(a* , —b*)), where q* = (a*,-b*), be a conjugate regular 
function. Its derivative at the point q* = (a*, —6*) is defined by the equation 

W ' Ag'^o Aq* 
CJV(A 9 *)#0 

provided this limit exists, when it is computed for nonsingular Aq* — (Aa*,— A6*). 
Analysis of this derivative leads to the equivalent complex pair representations 



dip 



(67) — = • '/ >> > 

(68) 



dq* \ da* ' da 

dip ( d4>2 d(pi 



dq* \ db* ' db* 

The equivalence of ( |67|) and (|68| ) presupposes the validity of the conjugate regular Cauchy- 
Riemann equations 



(69) 



(70) 



Q(Pl _ Q<P2 

da* db* 

d<f> 2 _ dcj>i 

da* ~ ~db~* 



Let us summarize our present findings concerning diffcrentiable functions. In Section 6.2 
we saw that in the complex plane a holomorphic or conjugate holomorphic function has 
an argument of the type a = (x, y) or a — (x, —y) when expressed as a real pair. In 
Chapter 4 and this section we have obtained holomorphic and conjugate holomorphic 
functions in the bicomplex space by applying the elementary functions of Chapter 3 to 
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the arguments p = (a, b), q = (a, b*), and q* = (a* , —b*). Looking at the form of these 
arguments it is hard to avoid the thought that there must be differentiable functions 
associated with each of the following arguments: 

(71) p=(a,b) , q=(a,b*) , r = (a*,b) , s = (a*,b*) 

(72) P u = (a,-b) , q" = (a,-b*) , r u = (a*,-b) , s u = (a*,-b*) 

Each argument of the second row is the bicomplex conjugate of the corresponding argu- 
ment of the first row. 

Bicomplex functions applied to the arguments (|7l|)-(|72]) are evidently of two main 
types: 

(73) m = {^XIH^X)) ,t=(a°,b«) 

(74) V(0 = (0i (a°, -&*), 02(^,-6°)) ,t u = (a D ,-&<>) 

The variable t should be replaced by p, q, r or s and either of the symbols ° and ° 
by a blank or a *. With the derivative defined as the limit of a fraction containing 
the appropriate argument the functions (|7^) are holomorphic if they fulfill complexified 
CR-equations of the form 

, , dfa = dfa 

V ' da a db° 

1 ' da? db° 



while the functions (|7J) are conjugate holomorphic if they fulfill 

(771 ^tl = d< ^ 2 

{ ' da° db" 

(7H) dfo = dfa_ 

K ' da a db" 

Equations ([75|)-([76|) correspond to the CR-equations (g), equations (|77|)-([78|) to the 
conjugate CR-equations (R). 



Example. Application of the bicomplex cosinus function (3.27) to r u = (a*,—b) yields 
cos r u = (cos a* ■ cosh b , sin a* ■ sinh b) . The component functions cos a* ■ cosh b and 
sin a* • sinhfr satisfy (j77|)~(|7g|) if a * is substituted for ° and a blank for °: 



oq>\ _ Q(P2 
da* db 

d(j>2 _ 34>i 
da* db 

□ 
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6.6 Bicomplex solutions of the Laplace equation 

A final remark about the solutions of the Laplace equation is in order. 

Let the function ip be given by ( |73| ) so that <j>\(a a ,b°) = tpi (x, y, z, u) + iip2(x, y,z,u) 
and 02(a D j b") = ips(x, V, z , u ) + iip4,(x, y, z, u). (The corresponding treatment of a func- 
tion of type (|74| ) is analogous.) If (f>\ and 4>2 are holomorphic or conjugate holomorphic 
in a and b, their component functions fulfill the two-dimensional Laplace equation: 

d*ipi+dfai = ,i= 1,2,3,4 
dl^i + dlipi = 

It follows that ^ fulfills the four-dimensional Laplace equation d^.ip+dyip+d^ip+d^ip = 0. 

Application of the differential operator to one of the Cauchy-Riemann equations 
j75|)-(|76|) and to the other yields after computing the sum or difference of the results 
either 

da? 2 db° 2 

or 

d 2 4>2 d 2 4>2 
da n2 db 02 

The complex pair components of all bicomplex holomorphic functions thus satisfy the 
two-dimensional, complexified Laplace equation. 

In fact, an alternative approach to the understanding of differentiable bicomplex 
functions is to take the previous two equations as starting point and combine them into 

(79) (^ + ^< O ) ^) = 

This is a generalization of equation (|l]). The operator to the left admits the factorization 

/ d 2 d 2 \ _ ( d_ d_ \ _ ( d_ _d_\ 

\da a2 + db° 2 ' ) ~ \da° ' db°) ® \da° ' ~ db° J 

which means that the solutions of 

(^•-^) e *-« - 

are solutions of (|7^) as well. Because these two equations are equivalent to (f75|)-(|76|) 
and (frjj)-(|78|), respectively, we see as in our treatment of the complex holomorphic 
and conjugate holomorphic functions in Section 6.1, that the corresponding bicomplex 
functions can be obtained by factorizing the appropriate version of the Laplace equation. 



= 
= 
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7 Conclusions 



In the course of this work we have frequently seen that a concept introduced or a result 
obtained has been the generalization of the corresponding concept or result in complex 
analysis. No doubt there are still plenty of such generalizations to be made, notably with 
regard to power series, Laurent series and analytic continuation. 

A somewhat separate matter is the development of a theory of tricomplex functions 
of the type fl = (\E , i,\E , 2), where \&i and ^2 are bicomplex functions of two bicomplex 
variables. In forming differentiable tricomplex functions one has at one's disposal all the 
holomorphic and conjugate holomorphic functions of Section 6.5. 

A corner stone of our work has been the application of complex pairs to represent 
quaternions and bicomplex numbers. With respect to the former essentially the same 
idea is found in the book jl9) by MacLane and Birkhoff, but instead of rendering a 
quaternion as a complex pair (a, b) they denote it by a + bj, which in our notation 
becomes (a, 0) + (6,0) x j. In |^| Price employs the notation a + i^b for the bicomplex 
number that we would write (a, 0)+jO(&, 0). All these notations seem rather cumbersome 
compared to the concise pair notation. 

In this context it may be noted that the practice of denoting complex numbers by 
pairs of real numbers is due to Hamilton JD|. His objective was to employ these algebraic 
couples or number-couples, as he called them, to remove the "metaphysical stumbling 
blocks" that he felt was connected with the imaginary unit y/—l. An account of Hamil- 
ton's thoughts on algebraic couples is found in |n| , the reading of which inspired the 
writing of this treatise. 
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Appendix 

The generalized Green theorem 

We shall give a proof of the generalized Green theorem as stated in Chapter 5. 

Theorem 5.4 Let <f>\ and (f> 2 be two complex functions of a and b such that they are 
holomorphic in a domain G C B. Also let S(T) be an orientable surface in G bounded 
by the closed curve T. Then 

(!) j^.*^.*.^*.*).*., 

Proof. Let the R 4 -representations of cpi and 4> 2 be 

(2) <Ai(a,&) = ipi(x,y,z,u) + iip 2 (x,y,z,u) 
(j) 2 (a,b) = ip 3 (x,y,z,u) +iip4:(x,y,z,u) 

a = x + iy , b = z + iu 

Because (f>\ and 4> 2 are holomorphic in G we get on application of (3.||) the following 
derivatives: 

(3) = d x ip 3 + id x ip4 

(4) ^ = d^+id.ifa 

ob 

Inserting (||) , (||) and (|J) into ([!]) , performing the complex multiplications and demanding 
the equality of real and imaginary parts change our proof obligation into: 
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(5) / (V'i dx — tp2 dy + ip3 dz — tp^ du) = 

[(d x 4> 3 - d z ipi) dxdz - (d x ips - d z ipi) dydu - 
(d x ip4 - d z ip 2 ) dxdu - (d x ip4 ~ d z ip 2 ) dydz)} 

( I i ) / ( < -j dx + tpi dy + -04 dz + ip 3 du) = 

[{d x ip 3 - d z ipi) dxdu + (d x ip 3 - d z ipi) dydz + 

S(T) 

(d x ip4 - d z ip 2 ) dxdz - (<9 X V>4 - d z ip 2 ) dydu)} 

The line and surface integrals of these formulas are of the type found at the right-hand 
sides of (5{[§) and (5|J|). 
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At this point we invoke the generalized Stokes theorem in R 4 g9| . It states that four 
real-valued, continuously differentiable functions A\, A 2 , A3, A4 of the variables x, y, z, u 
satisfy 

(7) J {Ax dx + A 2 dy + A3 dz + A A du) = 

I {A12 dxdy + A13 dxdz + A14 dxdu + A23 dydz + A24 dydu + A34 dzdu) 

JS(T) 

where S(T) is an orientablc surface in R 4 bounded by T and 

A12 = d x A 2 - d y Ai , A13 = d x A 3 - d z A x , A u = d x A 4 - d u A x 
A23 = d v A 3 - d z A 2 , A 24 = d y A 4 ~ d u A 2 , A 34 = d z A 4 - d u A 3 

Our intention is to show that (Q) reduces to (||) or (Q) when the A^-functions are chosen 
properly. Focusing on formula (0) first, we compare its left-hand side with the left-hand 
side of (0) and identify 

M = tpi , A 2 = -tp 2 , A 3 = i/j 3 , A 4 = -ij)4 

The functions <f>\ and $2 are holomorphic in a and 6, hence the ^-functions are contin- 
uously differentiable and satisfy the CR-equations (3j|) in the following way: 



d x ipi 


= dy1p2 


, d x ip2 = 


-Oylpl 


d z ipi 


= 9 u ip2 


, d z ip 2 = 


-d u ipi 




= 9 y ip4 


, 9 x tfj4 = 


~dy1p 3 


d z ip 3 


= d u i/; 4 


, d z ip4 = 


-d u if> 3 



We thus get 

A12 = -d x ij)2 - dytp! = 

Al3 = d x ip 3 - d z lp! 

A14 = -d x ip4, - d u 1p! = -8 x 1p4 + 8 z 1p2 

A 23 = dyip3 + d z il> 2 = ~(d x ip4- d z ip 2 ) 

A 2 4 = -d y ip4 + d u ip 2 = -d x ip3 + d z ipi 

A34 = -d z ip4 - d u ifj3 = 

These identities make the right-hand side of (0) equal to the right-hand side of (||), as 
required. 

Comparison of the left-hand side of (0) with the left-hand side of (^), in turn, tells 
us to choose 

Ai = ljj 2 , A 2 — tp\ , A3 — 1p4 , A4 = lf> 3 
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It is as straightforward as above to show that evaluated for this choice the expressions 
Aij make the right-hand side of (^) equal to the right-hand side of (^), which concludes 
the proof. 

□ 
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